BAB.ffi-PUB-12/018 
SLAC-PUB- 15181 



Measurement of 13{B — >■ Xs^), the B Xg^ photon energy spectrum, 
and the direct CP asymmetry in B ^ Xs-s^dl decays 

J. P. Lees, V. Poireau, and V. Tisserand 
Laboratoire d'Annecy-le-Vieux de Physique des Particules (LAPP), 
Universite de Savoie, CNRS/IN2P3, F-74941 Annecy-Le-Vieux, France 

J. Garra Tico and E. Granges 
Universitat de Barcelona, Facultat de Fisica, Departament ECM, E-08028 Barcelona, Spain 

A. Palano'^'' 

INFN Sezione di Bari"" ; Dipartimento di Fisica, Universitd di Bari'' , 1-70126 Bari, Italy 

G. Eigen and B. Stngn 
University of Bergen, Institute of Physics, N-5007 Bergen, Norway 

D. N. Brown, L. T. Kerth, Yu. G. Kolomensky, and G. Lynch 
Lawrence Berkeley National Laboratory and University of California, Berkeley, California 94720, USA 

H. Koch and T. Schroeder 
Ruhr Universitat Bochum, Institut fiir Experimentalphysik 1, D-44780 Bochum, Germany 

D. J. Asgeirsson, C. Hearty, T. S. Mattison, J. A. McKenna, and R. Y. So 
University of British Columbia, Vancouver, British Columbia, Canada V6T IZl 

A. Khan 

Brunei University, Uxbridge, Middlesex UBS 3PH, United Kingdom 

V. E. Blinov, A. R. Buzykaev, V. P. Druzhinin, V. B. Golubev, E. A. Kravchenko, A. P. Onuchm, 
S. L Serednyakov, Yu. L Skovpen, E. P. Solodov, K. Yu. Todyshev, and A. N. Yushkov 

Budker Institute of Nuclear Physics, Novosibirsk 630090, Russia 

M. Bondioli, D. Kirkby, A. J. Lankford, and M. Mandelkern 
University of California at Irvine, Irvine, California 92697, USA 

H. Atmacan, J. W. Gary, F. Liu, O. Long, and G. M. Vitug 
University of California at Riverside, Riverside, California 92521, USA 

C. Campagnari, T. M. Hong, D. Kovalskyi, J. D. Richman, and C. A. West 
University of California at Santa Barbara, Santa Barbara, California 93106, USA 

A. M. Eisner, J. Kroseberg, W. S. Lockman, A. J. Martinez, B. A. Schumm, A. Seiden, and L. Winstrom 
University of California at Santa Cruz, Institute for Particle Physics, Santa Cruz, California 95064, USA 

D. S. Chao, C. H. Cheng, B. Echenard, K. T. Flood, D. G. Hitlin, P. Ongmongkolkul, F. C. Porter, and A. Y. Rakitin 

California Institute of Technology, Pasadena, California 91125, USA 

R. Andreassen, Z. Huard, B. T. Meadows, M. D. Sokoloff, and L. Sun 
University of Cincinnati, Cincinnati, Ohio 45221, USA 

P. C. Bloom, W. T. Ford, A. Gaz, U. Nauenberg, J. G. Smith, and S. R. Wagner 
University of Colorado, Boulder, Colorado 80309, USA 



R. Ayacfl and W. H. Toki 



Colorado State University, Fort Collins, Colorado 80523, USA 

B. Spaan 

Technische Universitat Dortmund, Fakultat Physik, D-44S21 Dortmund, Germany 

K. R. Schubert and R. Schwierz 
Technische Universitat Dresden, Institut fiir Kern- und Teilchenphysik, D-01062 Dresden, Germany 

D. Bernard and M. Verderi 

Laboratoire Leprince-Ringuet, Ecole Polytechnique, CNRS/IN2P3, F-91128 Palaiseau, France 

P. J. Clark and S. Playfer 
University of Edinburgh, Edinburgh EH9 3JZ, United Kingdom 

D. Bettoni", C. Bozzi", R. Calabrese°^ G. Cibinetto"^ E. Fioravanti°^ 
I. Garzia°^ E. Luppi''^ M. Munerato°^ L. Piemontese" , and V. Santoro" 
INFN Sezione di Ferrara" ; Dipartimento di Fisica, Universitd di Ferrard' , F44IOO Ferrara, Raly 

R. Baldini-Ferroli, A. Calcaterra, R. de Sangro, G. Finocchiaro, 
P. Patteri, I. M. PeruzziQ M. Piccolo, M. Rama, and A. Zallo 
INFN Laboratori Nazionali di Frascati, F00044 Frascati, Raly 

R. Contri°^ E. Guido''^ M. Lo Vetere°^ M. R. Monge°^ S. Passaggio", C. Patrignani°^ and E. Robutti" 
INFN Sezione di Genova°^ ; Dipartimento di Fisica, Universitd di Genova'' , I-I6I46 Genova, Italy 

B. Bhuyan and V. Prasad 
Indian Institute of Technology Guwahati, Guwahati, Assam, 781 039, India 

C. L. Lee and M. Morii 
Harvard University, Cambridge, Massachusetts 02138, USA 

A. J. Edwards 
Harvey Mudd College, Claremont, California 91711, USA 

A. Adametz and U. Uwer 
Universitat Heidelberg, Physikalisches Institut, Philosophenweg 12, D-69120 Heidelberg, Germany 

H. M. Lacker and T. Lueck 
Humboldt-Universitdt zu Berlin, Institut fiir Physik, Newtonstr. 15, D-12489 Berlin, Germany 

P. D. Dauncey 
Imperial College London, London, SW7 2AZ, United Kingdom 

U. Mallik 

University of Iowa, Iowa City, Iowa 52242, USA 

C. Chen, J. Cochran, W. T. Meyer, S. Prell, and A. E. Rubin 

Iowa State University, Ames, Iowa 5001 1-3160, USA 

A. V. Gritsan and Z. J. Guo 
Johns Hopkins University, Baltimore, Maryland 21218, USA 

N. Arnaud, M. Davier, D. Derkach, G. Grosdidier, F. Le Diberder, A. M. Lutz, 
B. Malaescu, P. Roudeau, M. H. Schune, A. Stocchi, and G. Wormser 
Laboratoire de I'Accelerateur Lineaire, IN2P3/CNRS et Universite Paris-Sud 11, 
Centre Scientifique d'Orsay, B. P. 34, F-91898 Orsay Cedex, France 

D. J. Lange and D. M. Wright 



3 

Lawrence Livermore National Laboratory, Livermore, California 94550, USA 

C. A. Chavez, J. P. Coleman, J. R. Fry, E. Gabathuler, D. E. Hutchcroft, D. J. Payne, and C. Touramanis 
University of Liverpool, Liverpool L69 7ZE, United Kingdom 

A. J. Bevan, F. Di Lodovico, R. Sacco, and M. Sigamani 

Queen Mary, University of London, London, El iNS, United Kingdom 

G. Cowan 

University of London, Royal Holloway and Bedford New College, Egham, Surrey TW20 OEX, United Kingdom 

D. N. Brown and C. L. Davis 
University of Louisville, Louisville, Kentucky 40292, USA 

A. G. Denig, M. Fritsch, W. Gradl, K. Griessinger, A. Hafner, and E. Prencipe 
Johannes Gutenberg-Universitdt Mainz, Institut fiir Kernphysik, D-55099 Mainz, Germany 

R. J. BarlowH G. Jackson, and G. D. LafFerty 

University of Manchester, Manchester M13 9PL, United Kingdom 

E. Behn, R. Cenci, B. Hamilton, A. Jawahery, and D. A. Roberts 
University of Maryland, College Park, Maryland 20742, USA 

C. Dallapiccola 

University of Massachusetts, Amherst, Massachusetts 01003, USA 

R. Cowan, D. Dujmic, and G. ScioUa 

Massachusetts Institute of Technology, Laboratory for Nuclear Science, Cambridge, Massachusetts 02139, USA 

R. Cheaib, D. Lindemann, P. M. Pateli and S. H. Robertson 
McGill University, Montreal, Quebec, Canada H3A 2T8 

P. Biassoni'^^ N. Neri'', F. Palombo°^ and S. Stracka'^'' 
INFN Sezione di Milano"" ; Dipartimento di Fisica, Universitd di Milano^ , L20133 Milano, Italy 

L. Cremaldi, R. Godang0 R. Kroeger, P. Sonnek, and D. J. Summers 

University of Mississippi, University, Mississippi 38677, USA 

X. Nguyen, M. Simard, and P. Taras 
Universite de Montreal, Physique des Particules, Montreal, Quebec, Canada II3C 3J7 

G. De Nardo''^ D. Monorchio''^ G. Onorato"^ and C. Sciacca'^^ 
INFN Sezione di Napoli"" ; Dipartimento di Scienze Fisiche, 
Universitd di Napoli Federico if , 1-80126 Napoli, Italy 

M. Martinelli and G. Raven 
NIKHEF, National Institute for Nuclear Physics and High Energy Physics, NL-1009 DB Amsterdam, The Netherlands 

C. P. Jessop, K. Knoepfel, J. M. LoSecco, and W. F. Wang 
University of Notre Dame, Notre Dame, Indiana 46556, USA 

K. Honscheid and R. Kass 
Ohio State University, Columbus, Ohio 43210, USA 

J. Brau, R. Frey, M. Lu, N. B. Sinev, D. Strom, and E. Torrence 
University of Oregon, Eugene, Oregon 97403, USA 



E. Feltresi"^ N. Gagliardi"^ M. Margoni"^ M. Morandin", 



4 



M. Posocco", M. Rotondo", G. Siini", F. Si^lonetto°^ and R. Stroili"'' 

INFN Sezione dt Padova^ ; Diparttmento di Fisica, Universitd dt Padovci , 1-35131 Padova, Italy 

S. Akar, E. Ben-Haim, M. Bomben, G. R. Bonneaud, H. Briand, G. Calderini, 
J. Chauveau, O. Hamon, Ph. Leruste, G. Marchiori, J. Ocariz, and S. Sitt 
Laboratoire de Physique Nucleaire et de Hautes Energies, 
IN2P3/CNRS, Universite Pierre et Marie Curie-Paris6, 
Universite Dents Diderot-Parisl , F-75252 Paris, France 

M. Biasini''^ E. Manoni"^ S. Pacetti"^ and A. Rossi'''' 
INFN Sezione di Perugia'^ ; Dipartimento di Fisica, Universitd di Perugia'' , 1-06100 Perugia, Italy 

C. Angelini''^ G. Batignani''^ S. Bettarini''^ M. Carpinelli''''E!l G. Casarosa"^ A. Cervelli''^ F. Forti''^ 
M. A. Giorgi''^ A. Lusiani'^^ B. OberhoP^ E. Paoloni'^^ A. Perez'', G. Rizzo''^ and J. J. Walsh'' 
INFN Sezione di Pisa"^ ; Dipartimento di Fisica, Universitd di Pisa*" ; Scuola Normale Superiore di Pisa'^, 1-56127 Pisa, Italy 

D. Lopes Pegna, J. Olsen, A. J. S. Smith, and A. V. Telnov 

Princeton University, Princeton, New Jersey 08544, USA 

F. Anuhi", R. Faccini»^ F. Ferrarotto", F. Ferroni''^ M. Gaspero''^ L. Li Gioi'^, M. A. Mazzoni'^, and G. Piredda'' 

INFN Sezione di Roma"' ; Dipartimento di Fisica, 
Universitd di Roma La Sapienzd' , 1-00185 Roma, Italy 

C. Biinger, O. Griinberg, T. Hartmann, T. Leddig, H. Schroder^ C. Voss, and R. Waldi 
Universitdt Rostock, D-18051 Rostock, Germany 

T. Adye, E. O. Olaiya, and F. F. Wilson 
Rutherford Appleton Laboratory, Chilton, Didcot, Oxon, 0X11 OQX, United Kingdom 

S. Emery, G. Hamel de Monchenault, G. Vasseur, and Ch. Yeche 
CEA, Irfu, SPP, Centre de Saclay, F-91191 Gif-sur-Yvette, France 

D. Aston, D. J. Bard, R. Bartoldus, P. Bechtle, J. F. Benitez, C. Cartaro, M. R. Convery, J. Dorfan, 

G. P. Dubois-Felsmann, W. Dunwoodie, M. Ebert, R. C. Field, M. Franco Sevilla, B. G. Fulsom, A. M. Gabareen, 
M. T. Graham, P. Grenier, C. Hast, W. R. Innes, M. H. Kelsey, P. Kim, M. L. Kocian, D. W. G. S. Leith, P. Lewis, 

B. Lindquist, S. Luitz, V. Luth, H. L. Lynch, D. B. MacFarlane, D. R. MuUer, H. Neal, S. Nelson, M. Perl, 
T. Pulliam, B. N. Ratcliff, A. Roodman, A. A. Salnikov, R. H. Schindler, A. Snyder, D. Su, M. K. Sullivan, 
J. Va'vra, A. P. Wagner, W. J. Wisniewski, M. Wittgen, D. H. Wright, H. W. Wulsin, C. C. Young, and V. Ziegler 
SLAC National Accelerator Laboratory, Stanford, California 94309 USA 

W. Park, M. V. Purohit, R. M. White, and J. R. Wilson 

University of South Carolina, Columbia, South Carolina 29208, USA 

A. Randle-Conde and S. J. Sekula 

Southern Methodist University, Dallas, Texas 75215, USA 

M. Belhs, P. R. Burchat, and T. S. Miyashita 
Stanford University, Stanford, California 94305-4060, USA 

M. S. Alam and J. A. Ernst 
State University of New York, Albany, New York 12222, USA 

R. Gorodeisky, N. Guttman, D. R. Peimer, and A. Soffer 
Tel Aviv University, School of Physics and Astronomy, Tel Aviv, 69978, Israel 



P. Lund and S. M. Spanier 

University of Tennessee, Knoxville, Tennessee 37996, USA 



5 



J. L. Ritchie, A. M. Ruland, R. F. Schwitters, and B. C. Wray 
Umversity of Texas at Austin, Austin, Texas 78712, USA 

J. M. Izen and X. C. Lou 

University of Texas at Dallas, Richardson, Texas 75083, USA 

F. Bianchi"^ D. Gamba°^ and S. Zambito'''' 
INFN Sezione di Torino"" ; Dipartimento di Fisica Sperimentale, Universitd di Torino'' , T10125 Torino, Italy 

L. Lanceri"'' and L. Vitale"^ 
INFN Sezione di Trieste"' ; Dipartimento di Fisica, Universitd di Trieste'' , 1-34127 Trieste, Italy 

F. Martinez- Vidal and A. Oyanguren 
IFIC, Universttat de Valencia-CSIC, E-46071 Valencia, Spam 

H. Ahmed, J. Albert, Sw. Banerjee, F. U. Bernlochner, H. H. F. Choi, G. J. King, 
R. Kowalewski, M. J. Lewczuk, I. M. Nugent, J. M. Roney, R. J. Sobie, and N. Tasneem 
University of Victoria, Victoria, British Columbia, Canada V8W 3P6 

T. J. Gershon, P. F. Harrison, T. E. Latham, and E. M. T. Puccio 

Department of Physics, University of Warwick, Coventry CV4 7AL, United Kingdom 

H. R. Band, S. Dasu, Y. Pan, R. Prepost, and S. L. Wu 
University of Wisconsin, Madison, Wisconsin 53706, USA 
(Dated: July 23, 2012) 

The photon spectrum in B — >■ Xa'y decay, where Xs is any strange hadronic state, is studied using 
a data sample of (382.8±4.2) x 10® e+e" ~> r{AS) BB events collected by the BAB4R experiment 
at the PEP-II collider. The spectrum is used to measure the branching fraction B{B — > Xs'y) = 
(3.21 ± 0.15 ± 0.29 ± 0.08) x 10"* and the first, second, and third moments {E^) = 2.267 ± 0.019 ± 
0.032 ± 0.003 GeV, {{E^ ~ {E^)f) = 0.0484 ± 0.0053 ± 0.0077 ± 0.0005 GeV^ and {{E^ - {E^)f) = 
-0.0048 ± 0.0011 ± 0.0011 ± 0.0004 GeV^ for the range E^ > 1.8 GeV, where E-i is the photon 
energy in the B-meson rest frame. Results are also presented for narrower E^ ranges. In addition, 
the direct CP asymmetry Acp{B Xs+d^) is measured to be 0.057 ± 0.063. The spectrum itself 
is also unfolded to the B-meson rest frame; that is the frame in which theoretical predictions for its 
shape are made. 

PACS numbers: 13.25.Hw, 13.20.He, 12.15.Hh, ll.30.Er 



I. INTRODUCTION 



In the standard model (SM) the electromagnetic ra- 
diative decay of the b quark, b ^ or b dj , proceeds 
at leading order via the loop diagram shown in Fig. [T] 
resulting in a photon and a strange or down quark. The 
rate for 6 — >■ ^7 relative to 6 — >■ 57 is suppressed by 
a factor |Vtd/Vtsp where Vu and Vts are the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix elements. Interest 
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in these decays is motivated by the possibility that new 
heavy particles might enter into the loop at leading or- 
der, causing significant deviations from the predicted SM 
decay rates. There is an extensive theoretical literature 
evaluating the effects of new physics; some examples are 
given in references [IrS^]. New physics can also signifi- 
cantly enhance the direct CP asymmetry for b ^ and 
b^dj decay 

The hadronic processeses corresponding to the under- 
lying b and b ^ dj decays are B — >■ Xgj and 
B — )■ Xd'j- Here X^ and AT^ are any final state resulting 
from the hadronization of the sq or dq quark-level state, 
respectively, where q is the spectator from the B meson. 
These are predominantly resonances, including K*{892), 
K]^(1270) (A,) or p,uj (A^) and higher-mass states, but 
also non-resonant multi-hadron final states. Theoretical 
predictions for the rates of such exclusive decays suffer 
from large uncertainties associated with the form fac- 
tors of the mesons. In contrast, the inclusive hadronic 
rates T{B — )■ Xgj) and r{B — )■ Xd'y) can be equated 
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FIG. 1: The leading order Feynman diagram for the electro- 
magnetic radiative decay of the 6-quark in the SM. 



with the precisely-calculable partonic rates T(b — >■ 57) 
and T{b — >■ ^7) at the level of a few percent [ij] (quark- 
hadron duality), leading to significantly more accurate 
predictions. At next-to-next-to-leading order (up to four 
loops) the SM prediction for the branching fraction is 
B{B Xs-f) = 3.15 ±0.23 x IQ-'^iEj > 1.6 GeV) [H]. 
Measurements of the inclusive rates and asymmetries are 
therefore powerful probes of physics beyond the standard 
model. 

The shape of the photon energy spectrum is deter- 
mined by the strong interaction of the b quark within the 
B meson and by the hadronization process. The Fermi 
motion of the quark within the B meson and gluon ra- 
diation lead to an E-y distribution, in the _B-mcson rest 
frame, that is peaked in the range 2.2 to 2.5 GeV, with 
a kinematic limit at 7715/2 ss 2.64 GeV and a rapidly- 
falling low-e nerg y tail. The shape is insensitive to non- 
SM physics [l6l [iTj . and can therefore provide infor- 
mation about the strong-interaction dynamics of the b 
quark. Heavy Quark Effective Theory (HQET) [H, [13- 
I22 ] has been used most extensively to describe these dy- 
namics. The shape of the photon spectrum provides in- 
formation on parameters of this theory related to the 
mass and momentum of the b quark within the B meson; 
the definitions and hence the values of these parameters 
differ slightly between the "kinetic scheme" [23] and the 
"shape function scheme" 2J|. The Heavy Flavor Aver- 
aging Group (HFAG) [25| has computed world-average 
values of the parameters in the kinetic scheme based on 
previous measurements of the inclusive semileptonic B- 
meson decay B — ^ Xcd^ {£ = e 01 ji) and of -B — t- ^^7. 
HFAG has also translated those values to the shape func- 
tion scheme. These parameters can be used to reduce 
the error in the extraction of the CKM matrix elements 
\Vcb\ and \Vub\ from the inclusive semileptonic decays, 
B XJi' and B XJf [IB-Ill. The B -J> ^^7 spec- 
tral shape may also be compared to predictions in the 
framework of Dressed Gluon Exponentiation [s^ . 

The inclusive decay B — >■ ^^7 was first measured by 
the CLEG collaboration p3l[j33l an d has been subse- 
quently studied by the ALEPH lH, Belle and 
BaBAR [41-43] collaborations. All measurements have 
been made with B mesons produced in e^e~ collisions. 
The theoretical predictions, which assume that the mea- 



surement is inclusive so that quark-hadron duality holds, 
are made in the _B-meson rest frame for photons with 
-E-y > 1.6 GeV. This means that ideally the measure- 
ment is made for all X^ final states and for all photons 
> 1.6 GeV. The experimental challenge is to make 
the measurement as inclusive as possible while suppress- 
ing backgrounds from other processes producing photons 
or fake photons. The backgrounds arise from continuum 
events {e'^e^ to qq or r+T~ pairs, where q = u, d, s 
or c), with the photon coming from either a tt*^ or 77 de- 
cay or from initial state radiation, and from other BB 
processes. The BB background arises predominantly 
from 7r° or 77 decay but also from decays of other light 
mesons, misreconstructed electrons and hadrons. It is 
strongly dependent on photon energy and rises steeply 
at lower Ej. This places a practical lower limit for Ej 
on the experimental measurements; measurements have 
been made to date with E-y > 1.7, 1.8, and 1.9 GeV. 

Three experimental techniques have been pursued. 
They differ in the extent to which the final state is re- 
constructed. The first is the fully inclusive technique in 
which neither the X^ from the signal B nor the recoil- 
ing B meson is reconstructed. (Charge conjugates are 
implied throughout this paper.) The second is the semi- 
inclusive technique, in which as many exclusive X^ final 
states as possible are reconstructed and combined. The 
recoiling B meson is not reconstructed. The third is the 
reconstructed recoil-i? technique, in which inclusive B 
events are tagged by fully reconstructing the recoiling 
B mesons in as many final states as possible, but X^ is 
not reconstructed. Each of the techniques has different 
strengths and weaknesses. 

If the X^ is not reconstructed, the sample includes all 
X^ final states, but there are significant backgrounds from 

other BB decays that must be estimated. It also includes 
Xj^ states from the Cabbibo-suppressed b dj process. 
These can be subtracted by assuming the b dj pho- 
ton spectrum to have a similar shape to the b ^ 
photon spectrum, but scaled by the ratio of the CKM 
elements {\Vtd\ / \Vts\f = 0.044 ± 0.003. This is believed 
to be a valid assumption. Also, if the X^ is not recon- 
structed then the signal B cannot be reconstructed. The 
B mesons have a small momentum in the T(45') rest 
frame. As the B meson is not reconstructed, the direction 
of the momentum is not known. This leads to a Doppler 
smearing of the photon energy. This effect, along with 
the detector resolution, must be corrected for or unfolded 
in order to compare to predictions made in the S-meson 
rest frame. Quantities measured in the T{AS) rest frame 
(CM frame) such as the photon energy E* are denoted 
with an asterisk. 

No semi-inclusive measurement to date has recon- 
structed more than about 60% of X^ decays, due to the 
high combinatoric background for higher multiplicity de- 
cays. Uncertainties in modeling the mix of X^ final states 
result in significant efficiency uncertainties, as well as a 
large uncertainty in correcting for the final states that are 
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not reconstructed. However, the reconstruction of the 
implies that the signal B can be fully reconstructed, pro- 
viding kinematic constraints to strongly suppress back- 
grounds, allowing the measurement to be made directly 
in the _B-meson rest frame. 

In the reconstructed recoil- i? technique, only about 1% 
of B's can be fully reconstructed, due to the presence of 
neutrinos in semileptonic decays and combinatoric back- 
grounds to higher multiplicity decays. This severely lim- 
its the statistical precision, but does allow the measure- 
ment to be made in the _B-meson rest frame. 

This paper reports a fully inclusive analysis that su- 
persedes the previous BaBAR fully- inclusive result [4^ . 
which is based on a smaller data sample. The E* pho- 
ton spectrum is measured in B — ^ Xg+d^J decays. It is 
used to measure the branching fraction B{B — > Xg'y) 
for > 1.8 GeV and for narrower energy ranges. The 
effects of detector resolution and Doppler smearing are 
unfolded to provide an E^ photon spectrum in the B- 
meson rest frame that can be used to fit to theoretical 
predictions for the spectral shape. The unfolded spec- 
trum is also used to measure the first, second and third 
moments, given respectively by 

El = (E^) 

E, = {iE,~{E,))^) (1) 
Es = {{E^ - {E^)f) . 

Although the SM predicts quite different asymmetries 
for B — >• Xgj and B — )• Xdj, the X^ and X^ final states 
cannot be distinguished in the fully inclusive technique. 
Hence the B — > X^'y contribution to the fully-inclusive 
measurement cannot be corrected for, and only the com- 
bination Acp {B — > Xs+dl) can be measured: 

r(6 ^ S7 + 6 ^ ^7) - r(5 ^ S7 + 6 ^ d7) 
Acp = = = = — . 

r(6 ^ S7 + 6 ^ ^7) + r(6 ^ S7 + 6 ^ ^7) 

This asymmetry is approximately 10^® in the SM, with 
nearly exact cancellation of opposite asymmetries for 
6 — )■ 57 and b — d-f. Acp{B Xg+dl) is sensitive to dif- 
ferent new physics scenarios than Acp {B — )■ X^j) pT| . 
Thus measurements of t his j oint asy mmetry complement 
those of Acp in 6 — 57 |32l . [36l . |44 145| to constrain new 
physics models. 

II. DATASETS, DETECTOR, SIMULATION 
AND SIGNAL MODELS 

The results presented are based on data samples of 
e+e" — >■ T{AS) — >■ BB collisions collected with the 
BaBAR detector at the PEP-II asymmetric-energy e~^e~ 
collider. The on-resonance integrated luminosity is 
347.1 fb~^, corresponding to 382.8 million BB events. 
The continuum background is estimated with an off- 
resonance data sample of 36.4 fb"^ collected 40MeV be- 
low the T{AS) resonance energy. 



The BaBAR detector is described in detail in refer- 
ence [46J. Charged-particle momenta are measured with 
a 5-layer, double-sided silicon vertex tracker (SVT) and 
a 40- layer drift chamber (DCH) inside a 1.5-T super- 
conducting solenoidal magnet. A high resolution total- 
absorption electromagnetic calorimeter (EMC), consist- 
ing of 6580 CsI(Tl) crystals, is used to measure local- 
ized electromagnetic energy deposits and hence to iden- 
tify photons and electrons. The EMC energy resolution 
for high-energy photons in the current measurement is 
about 2.6%. A ring-imaging Cherenkov radiation detec- 
tor (DIRC), aided by measurements of ionization energy 
loss, dE/dx, in the SVT and DCH, is used for parti- 
cle identification (PID) of charged particles. Muons are 
identified in the instrumented flux return (IFR), which 
consists of 18 layers of steel interleaved with single-gap 
resistive-plate chambers. For the last 38% of the data col- 
lected, 1/3 of these chambers in the central region of the 
detector were replaced by 12 layers of limited-streamer 
tubes, interspersed with 6 layers of brass (to increase ab- 
sorption) . 

The BaBAR Monte Carlo (MC) simulation, based on 
GEANT4 47], EVTGEN Q and JETSET [H, is 
used to generate samples of B~^B~ and B^B^ , qq (where 
g is a M, d, s or c quark), r+r^, and signal events (BB 
events in which at least one B decays to ^^7). To model 
beam backgrounds, each simulated event is overlaid with 
one of a set of of random background data events col- 
lected using a periodic trigger. 

The signal models used to determine selection efficien- 
cies are based on QCD calculations of references [l^ (ki- 
netic scheme) and [13] (shape function scheme) and on 
an earlier calculation by Kagan and Neubert Q ("KN"). 
Each model uses an ansatz for the shape that is con- 
strained by calculations of the first and second moments 
of the spectra. The models approximate the hadronic 
mass (toXs ) spectrum, which contains a number of over- 
lapping resonances, as a smooth distribution. This is rea- 
sonable, except at the lowest masses, where the i4'*(892) 
dominates the spectrum. Hence the portion of the mx, 
spectrum below 1.1 GeV/c^ is replaced by a Breit-Wigner 
if* (892) distribution, normalized to yield the same frac- 
tion of the integrated spectrum. A particular signal 
model is defined as the theoretical spectrum for specific 
HQET parameters, with this K*{892) at low nix,- The 
photon energy in the _B-meson rest frame is related to 
mx, via 

E 

"^x =m|-2m_B^. (2) 

High-statistics MC signal samples for the non-i4r*(892) 
part of the spectrum are generated uniformly in E-y- sep- 
arately for each of the two B-meson charge states - and 
then weighted according to any particular model of in- 
terest. 

Monte Carlo samples of B^B^ and B^B'^ events are 
needed for background evaluation. They are produced, 
with nearly three times the effective luminosity of the 



8 



data sample, and include all known B decays, except for 
events in which either B decays via B — >■ Xg+dl- Monte 
Carlo samples of continuum events {qq, separately for cc 
and for the light quarks, and t~^t^) are used to optimize 
the event selection criteria, but are not otherwise relied 
upon. 



III. ANALYSIS OVERVIEW 

The event selection is described in detail in Sec. IIVI 
The analysis begins by selecting hadronic events. A 
high-energy photon, characteristic of B — > Xsj decays, 
is then required, while photons from tt" and r] decays 
are vetoed, reducing both the continuum and BB back- 
grounds. The background from continuum events is sig- 
nificantly suppressed by charged lepton tagging (requir- 
ing a high-momentum lepton, as would be expected from 
the semileptonic decay of a i? meson) and by exploit- 
ing the more jet-like topology of the qq or t+t~ events 
compared to the isotropic BB decays. 

The continuum MC simulation does not adequately 
model the actual continuum background, primarily be- 
cause it omits QED and two-photon processes. Hence the 
continuum background is estimated with off-resonance 
data (Sec.|V]), which limits the statistical precision of the 
signal-yield measurement. However, the continuum sim- 
ulation is used to optimize some of the event-selection 
criteria (which must be done without reference to ac- 
tual data). After preliminary event selection, which re- 
duces the unmodeled backgrounds, a simple scaling of the 
continuum-MC predictions adequately models the event- 
yield distributions relevant for optimization. 

The lepton tagging and event topology criteria do not 
substantially reduce the BB background relative to the 
signal, as these processes have similar characteristics. 
The remaining BB background is estimated using MC 
simulation. There are several different i?-meson decays 
that contribute. Section I VII describes how each signifi- 
cant component is compared to an independent data con- 
trol sample and weighted to replicate those data. The un- 
certainty in these weighting procedures is the dominant 
source of systematic uncertainty. 

After the event selection, the continuum and 
reweighted BB backgrounds are subtracted from the 
on-resonance data sample, resulting in the raw B — > 
Xs+dl photon spectrum (Sec. IVHp . The analysis was 
done "blind" in the range of reconstructed photon en- 
ergy E* from 1.8 to 2.9 GeV; that is, the data were not 
looked at until all selection requirements were set and the 
corrected backgrounds determined. The choice of signal 
range is limited by high BB backgrounds at low E* . The 
regions 1.53 < E* < 1.8 GeV and 2.9 < E* < 3.5 GeV 
are dominated hy BB and continuum backgrounds re- 
spectively. They provide control regions to validate the 
background estimation for the signal region. 

The raw spectrum is used to extract the direct CP 
asymmetry (Sec. IIX[) and the partial branching fraction 



for 1.8 < E^ < 2.8 GeV (Sec. El. Finally, in Sec.|XT]the 

effects of detector resolution and Doppler smearing are 
unfolded in order to measure the shape of the photon 
energy spectrum in the _B-meson rest frame. 

IV. EVENT SELECTION 

The event selection was developed using MC samples 
of signal and background events. The model used for 
signal simulation, as defined in Sec. [TTl is based on a KN 
spectrum with rrib = 4.65GeV/c^. 

A. Selection of Hadronic Events 

For each event, the analysis considers good-quality re- 
constructed tracks, which have momenta transverse to 
the beam direction of at least O.lGeV/c and originate 
from the vicinity of the interaction point (point of closest 
approach within 10 cm along the beam axis and 1.5 cm 
in the transverse plane), and EMC clusters of at least 
30 MeV in the laboratory frame. Hadronic events are se- 
lected by requiring at least three reconstructed charged 
particles and the normalized second Fox- Wolfram mo- 
ment [131 -^2 to be less than 0.90. To reduce radia- 
tive Bhabha and two-photon backgrounds, the number 
of charged particles plus half the number of photons with 
laboratory-frame energy above 0.08 GeV is required to be 
greater than 4.5. 

B. Requirements on the High-Energy Photon 

The photon selection requires at least one photon can- 
didate with 1.53 < E* < 3.5 GeV in the event. A photon 
candidate is a neutral EMC energy cluster with a lateral 
moment consistent with that of a single photon [5l| . The 
latter requirement rejects most background from neutral 
hadrons, which at these energies is dominated by antineu- 
trons that annihilate in the EMC. The photon location 
is assigned at a depth of 12.5 cm in the EMC, where it 
is required to be isolated by 25 cm from any other en- 
ergy deposit (the lateral dimensions of the crystals are 
approximately 5 cm by 5 cm) . The cluster must also be 
well-contained in the calorimeter (—0.74 < cos9^ < 0.94, 
where 6^ is the laboratory-frame polar angle with respect 
to the direction of the electron beam) . A likelihood vari- 
able (L^a) based on the energy profile of the EMC cluster 
is used to suppress the contribution of tt^'s in which the 
two daughter photons are not resolved. The requirement 
on L^rO retains essentially all isolated high-energy pho- 
tons. These photon quality criteria are determined from 
studies of photons in fi^j events and of p's (p's) from A 
(A) decays. (Antiprotons are used to estimate the detec- 
tor response to background antineutrons.) 

High-energy photons that are consistent with originat- 
ing from tt" — 77 or 77 — 77 decays are vetoed if the 
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other 7r° or r] daughter is found. For the tt°{i]) veto, com- 
binations are formed of the high-energy photon with all 
other photon candidates that have laboratory-frame en- 
ergy greater than 30 (230) MeV; it is required that the 
invariant mass not lie within a window around the nom- 
inal TT^irj) mass, 115(508) < m^^ < 155(588) MeV/c^. 

The simulated distributions of signal and background 
at this stage of the event selection are shown in Fig. ^a) . 
The cumulative signal efficiency up to this point is ap- 
proximately 50%, while 1.6% of continuum and 0.4% of 
BB backgrounds are retained. The remaining contin- 
uum background arises predominantly from unvetoed 7r° 
and T] decays, or initial state radiation in qq events. The 
BB background is also dominated by unvetoed decays of 
7r°(r/) from B — X7r°(?7), but also has a significant con- 
tribution from misidentified electrons, and smaller com- 
ponents from antineutrons and radiative uj and r/ decays. 

C. Lepton Tagging 

About 20% of B mesons decay semileptonically to ei- 
ther an electron or muon, predominantly via B — >■ Xcii^. 
An additional 4% of B decays result in an electron or 
muon via B — > XcTi^. Since the tagging lepton comes 
from the recoiling B meson, this requirement does not 
compromise the inclusiveness of the B — >■ Xg^ selection. 

Electrons are identified with a likelihood algorithm 
that incorporates properties of the deposited energy and 
shower shapes of the EMC clusters, the Cherenkov angles 
associated with the charged particle passing through the 
DIRC, and the dE/dx energy loss of the track. Muons 
are identified using a neural-network selector contain- 
ing variables that discriminate between muons and elec- 
trons, primarily through differences in EMC energy depo- 
sition, and those which discriminate between muons and 
hadrons, mainly through differences in IFR signatures. 

The left plots of Fig.[3]show that leptons from hadronic 
decays in continuum events tend to be at lower momen- 
tum. Hence the tagging lepton is required to have mo- 
mentum pI ^ > 1.05GeV/c. As seen in the right plots 
of Fig. [31 additionally requiring the cosine of the CM- 
frame angle between the lepton and the high-energy pho- 
ton cos 0*^ > —0.7 removes more continuum background, 
in which the lepton and photon candidates tend to be 
back-to-back. The peak at cos6'*^ ~ 1.0 for electrons in 
continuum events arises predominantly from 7r°(77) — ^ 77 
decays in which one photon satisfies the high-energy pho- 
ton requirements and the other converts to an e~^e~ pair. 
The peaks at cosO*^ ~ —1.0 for the BB background 
arise from B decays in which the photon and lepton come 
from the same B. A similar smaller peak for muon tags 
in signal events is due to pions faking the muon signa- 
ture. These tag selection requirements are designed as 
a loose preselection; a more stringent tag discrimina- 
tion is achieved by the multivariate selectors described 
in Sec. HVP] 

The presence of a relatively high-energy neutrino in 



semileptonic B decays is exploited by requiring the miss- 
ing energy of the event {E^:^^J to be greater than 0.7 GeV. 
The lepton-tag requirements retain approximately 12% 
of signal and BB background events after the photon 
selection, while retaining only 2.2% of continuum back- 
grounds. 

D. Event Topology Requirement 

As the continuum backgrounds are different for elec- 
tron and muon tags, each sample is divided according to 
the tag. For each lepton type the continuum backgrounds 
are then further suppressed by combining the p* ^ and 
cos 9*g for the leptons with event topology variables into 
a neural network (NN) discriminant. 

Several alternative choices of input variables were con- 
sidered. For each alternative, the electron and muon 
NN's are trained, and the requirements on their output 
parameter optimized (see below) . The choice of variables 
is designed to minimize the total error on the branching 
fraction and spectral moment measurements, based on 
combining in quadrature preliminary estimates of statis- 
tical, systematic and model-dependence errors. The lat- 
ter refers to a variation of the event-selection efficiency 
with the choice of MC spectrum ("model" in the sense of 
Sec. |TT]) used to compute it. It arises primarily from the 
increase in efficiency as a function of E* ; the stronger this 
trend, the larger the model-dependence uncertainty. The 
selection strategy aims for best signal precision, while 
minimizing the dependence of efficiency on E*. Since 
the backgrounds rise sharply as E* decreases, it is im- 
possible to completely eliminate the E* dependence. Of 
several multivariate discriminants (with different sets of 
input variables) which were found to give approximately 
the same signal precision, the one resulting in the least 
E* dependence was chosen. 

The eight topology variables chosen for the NN include 
R2/R2, where i?2 is the normalized second Fox- Wolfram 
moment calculated in the frame recoiling against the pho- 
ton, which for ISR events is the qq rest frame. Also 
included are three momentum-weighted polar-angle mo- 
ments, Lj/LQ,j — 1,2,3, where 

Lj = ^ |pi|| cos6'i|-' . (3) 

i 

Here pj and cos 9i are the momentum and angle, respec- 
tively, of the ith reconstructed particle with respect to 
the high-energy photon axis in the recoil frame. Sum- 
mation over i includes every reconstructed charged and 
neutral particle except the high-energy photon. The last 
four topology variables are derived from the eigenvalues 
(Ai , A2 , A3) and eigenvectors of the momentum tensor [s^ 

Tjnva _ EjP"pr/ip*i 

where is the nth component of the ith reconstructed 
particle's 3-momentum in the recoil frame. The high- 
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FIG. 2: Estimated signal and background yields vs. photon energy in the CM frame based on MC simulation, at two stages of 
the event selection: (a) after requiring an unvetoed high-energy photon (logarithmic scale); (b) after all selection requirements 
(linear scale) . The three contributions are shown cumulatively. The signal distribution is for a KN model with mi, — 4.65 GeV/ , 
while the continuum distribution has been scaled as described in Sec. IIIII 




CM Muon Momentum (GeV/c) cos 6* (|i,y) 

FIG. 3: Lepton distributions from MC simulation, after the photon selection requirements but before applying lepton tag and 
NN criteria. Plots (a) and (b) are for electron tags, plots (c) and (d) for muon tags. Plots (a) and (c) show the CM-frame 
momentum distributions, with vertical lines indicating the minimum selection requirements. Plots (b) and (d) show the cosine 
of the CM angle between the lepton and the high-energy photon, after applying the momentum criteria; the vertical lines show 
the minimum requirement on this quantity. The signal (black dots) is from a KN model with mj, = 4.65 GeV/ c^. The BB 
background (solid blue histogram) and continuum background (dashed red) are from the MC simulations. Each distribution is 
separately normalized to best illustrate its behavior. 
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electron-tagged events and (b) muon-tagged events. The vertical lines show the minimum requirement on this quantity. The 
signal is from a KN model with rrib — 4.65GeV/c^. The continuum is from the MC simulation. Normalizations are arbitrary. 



energy photon candidate is excluded. The derived quan- 
tities used as NN inputs are: 

A" max(Ai, A2, A3) 

A^"^ = A1A2 + A2A3 + A3A1 

A'^'* = A1A2A3 

V^"^ — z-component of Vmax , 

where Vmax is the eigenvector associated with the largest 
eigenvalue and z is the electron beam direction. 

The electron and muon NN's are trained with MC 
samples of continuum and signal (KN model with = 
4.65GeV/c^) events that contain a photon with energy 
in the range 1.9 < E* < 2.7 GeV. The BB background 
simulation sample is excluded from the training because 
this sample is used for background subtraction and is 
topologically very similar to the signal. Training with 
background and signal samples normalized to the ex- 
pected event yields at this stage of the event selection 
provides slightly better statistical precision for signal (see 
Sec. lIVEp than does training with background and sig- 
nal samples with the same normalization. For an NN 
with equally-normalized training samples, the NN output 
distributions would peak toward and 1, respectively, 
for background- and signal-like events. Neural-network 
training based on expected event yields, however, pro- 
duces output distributions that are qualitatively differ- 
ent, as demonstrated in Fig. 21 which shows the out- 
put distributions for signal and continuum events, sepa- 
rated according to lepton tag. Events with an electron 
(muon) tag are required to have an NN output greater 
than 0.53 (0.47). This selection accepts 42% of signal 
events (1.8 < E* < 2.8 GeV) that have passed the photon 
and lepton selection requirements while retaining 1.7% 
of continuum and 27% of BB background. Events with 
more than one photon candidate after the NN require- 
ment are discarded (0.16% of signal events). 



E. Optimization of the Event Selection 

The optimization for the selection criteria was per- 
formed iteratively on five variables: the two NN outputs 
fSec. lIV D)) . the minimum energy of the lower energy pho- 
ton in the 7r° and rj vetoes (Sec. lIVBp and the missing 
energy (Sec. lIVCp . The figure of merit (FOM) is the 
anticipated ratio of the signal yield to its statistical un- 
certainty for E* between 1.8 and 2.8 GeV, taking into 
account the limited size of the off-resonance sample used 
for continuum subtraction: 

S 

Statistical FOM = , . (5) 

v/S + B + (C//off ) 

Here S, B and C are the estimated yields in the on- 
resonance data of signal, BB background and contin- 
uum background events, respectively (after event selec- 
tion) , based on MC simulation, and /off is the fraction of 
total luminosity accumulated off-resonance, £off/('Con + 
Cog) = 0.0949. 

The selection criterion for each of the five variables was 
optimized in turn, while holding the criteria for the others 
fixed, and the process repeated until a stable optimal 
selection was found. 



F. Overall Signal Efflciency 

The probability that a signal event is observed and sur- 
vives the event-selection process is approximately 2.5%, 
while only 0.0005% of the continuum and 0.013% of the 
BB backgrounds remain in the sample. Figure [2jb) 
shows the expected signal and background distributions 
after all selection criteria. 

The photon spectrum is measured in bins of recon- 
structed E*. Hence the signal efficiency is presented 
here in terms of that quantity. The selection efficiency 
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FIG. 5: Combined acceptance and event-selection efficiency 
vs. measured E* for a KN model with mt = 4.65GeV/c^. 
Uncertainties are from MC statistics. Corresponding efficien- 
cies for a kinetic-scheme model with parameters set to HFAG 
world-average values are within 1% (relative) of the values 
plotted. 



for MC signal events, i.e., the fraction of the events in 
a given range of E* that survive all the selection cri- 
teria described above, is calculated in 100-MeV bins of 
reconstructed E* and also for wide ranges (such as 1.8 
to 2.8 GeV). The overall signal efficiency also includes an 
acceptance component, the probability for the photon to 
enter the fiducial region of the EMC. This is available 
only as a function of true E* (the photon energy before 
resolution smearing), since reconstructed E* is defined 
only for accepted photons. However, because the varia- 
tion of the acceptance efficiency is weak, it can be com- 
bined with the selection efficiency to provide an overall 
efficiency in bins or ranges of reconstructed E*. Figure [S] 
shows the result. 



CONTINUUM BACKGROUNDS 



The continuum background is estimated using off- 
resonance data scaled according to the ratio of the lu- 
minosity of the on- and off-resonance data sets. Since 
continuum data are collected 40MeV below the T(4S') 
resonance, the center-of-mass energy is 0.4% lower than 
the center-of-mass energies for a typical BB event. In 
order to account for this difference, the energy of a high- 
energy photon candidate in off-resonance data is scaled 
by nir{4S)/y/SoSi where mr(4S) and ^Soff are the mass 
of the T{4S) system and the center-of-mass energy of the 
ofF-resonance data event, respectively. 



VI. BB BACKGROUNDS 
A. Overview 

The background from non-signal BB events arises ei- 
ther from real photons from the decays of low-mass 
mesons (with 7r° and rj responsible for most of the back- 
ground) or from other particles faking photons. 

The BB background remaining after event selection is 
estimated using the MC simulation as an approximate 
starting point. Various control samples are then used to 
correct most of the significant components of this back- 
ground according to data/MC yield ratios measured as 
a function of appropriate kinematic variables. The cor- 
rections are applied in 100-MeV bins of E*. The uncer- 
tainties of these factors (along with small uncertainties 
from MC statistics) constitute the BB systematic errors. 
These can be highly correlated between E* bins. The re- 
mainder of Sec. IVII details the individual corrections, as 
well as a more global correction to the lepton-tagging 
efficiency. 

The event simulation tells us the true (generated) par- 
ticle that most closely corresponds to the reconstructed 
high-energy photon candidate. This allows the catego- 
rization of selected events according to the origin of that 
candidate. Table H] lists the MC fractions by category 
and the corresponding correction factors averaged over 
two broad E* intervals, covering the BB control region 
and the signal region. 



B. tt" and rj Corrections 

About 80% of MC-predicted B'B background in the 
signal region arises from B — ^ Xtt'^It]) with Tr^(j]) — > 77. 
This contribution is dominated by highly-asymmetric 
7r'^(r/) decays, in which a second photon has much lower 
energy than the selected high-energy photon. To cor- 
rect MC predictions for these inclusive B decays in the 
phase space region selected for the B ^ X^^ analysis, 
inclusive tt^ and 77 samples are defined by applying the 
same selection criteria but omitting the tt" and 77 vetoes. 
To enhance statistics for these studies the minimum re- 
quirement on E* is relaxed from 1.53 GeV to 1.03 GeV, 
and for ry's the minimum laboratory- frame energy for the 
low-energy photon is relaxed from 230 MeV to 75 MeV. 



1. Scaling of MC and rj yields to data 

The yields of tt'^{'I]) are measured in bins of -E*o(^) by 

fitting the distributions of 77 mass (m) in simulated BB 
background, on-resonance data and off-resonance data. 

The signal shape for tt'^ is the sum of two Gaussian 
functions (Gi and G2) with different means {fj,i and ^2) 
and rms widths (cti and (72) plus a low-mass power-law 
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TABLE I: The BB background composition after all selection cuts, according to the BABAR Monte Carlo simulation, and 
the correction factors determined for each component. Classification is according to the true MC particle associated with the 
high-energy photon, and to the true parent of that particle. The "B" category under "Parent" corresponds to high-energy 
photons from final-state radiation. The "Other" category consists of hadrons other than n's. The "None" category consists of 
backgrounds unassociated with the primary event, mostly from out-of-time Bhabha-scattering events; such "photons" appear 
in the simulation via the beam-background mixing described in Sec. [Til While all numbers are actually computed and applied 
in 100-MeV bins of E*, they are illustrated here for the overall signal region (1.80 to 2.80 GeV) and BB control region (1.53 
to 1.80 GeV). The "Subsection" column refers to where each correction is discussed. Note that the vr" and rj correction factors 
implicitly include the tagging efficiency correction described in Sec. IVIHI not included elsewhere in this table. 



MC Category 


1.53 to 1.8 GeV 


1.8 to 2.8 GcV 


Subsection 


Particle 


Parent 


MC Fraction 


Corr. Factor 


MC Fraction 


Corr. Factor 


Photon 




0.5390 


1.05 


0.6127 


1.09 


IVIBI 




V 


0.2062 


0.79 


0.1919 


0.75 


IVIBl 






0.0386 


0.80 


0.0270 


0.80 


|VIC| 




ri' 


0.0112 


0.52 


0.0082 


1.13 


|VIC| 




B 


0.0362 


1.00 


0.0194 


1.00 


IVIFl 






0.0061 


1.00 


0.0071 


1.00 


IVIHl 




e± 


0.0967 


1.07 


0.0619 


1.07 


IVIDl 




Other 


0.0035 


1.00 


0.0032 


1.00 


IVIHl 




Total 


0.9375 




0.9315 






e± 


Any 


0.0411 


1.65 


0.0333 


1.68 


IVIDl 


n 


Any 


0.0170 


0.35 


0.0243 


0.15 


jVIEj 


Other 


Any 


0.0029 


1.00 


0.0028 


1.00 


IVIHl 


None 




0.0015 


1.00 


0.0079 


1.00 





tail (parameters p and A) : 

A [/iGi(m) + (1 - /i)G2(m)] m > mo 



f{m) = 



B 



mo 



m < Too , 



(6) 

where too = {pi — Acti), A and /i govern the normal- 
izations of the two Gaussian functions, and B is set by 
requiring continuity at m = too. The signal shape for rj 
is a single Gaussian with two such power-law tails with 
separate parameters. 

The fit is carried out in several stages. First, a sig- 
nal shape is determined for BB MC events in which the 
reconstructed 77 pair derive from a true 7r° or 77. For 
purposes of this study, these events are termed "signal" . 
Next, for both MC and on-resonance data events, the 
mass spectrum of all 77 pairs which include the high- 
energy photon is fit in the Tr^ii]) mass region to signal 
plus a background shape, with some signal tail parame- 
ters fixed to their values from the signal-only fit. This 
procedure is validated by comparing the extracted signal 
yields from this MC fit to those of true signal: averaging 
the absolute values of the differences over energy bins, the 
agreement is 1.3% of the yield for tt" and 2.1% for 77. The 
fits to on-resonance data are shown in Figs. [6] and [71 Fi- 
nally, the off-resonance data are fit with all signal-shape 
parameters fixed to their on-resonance fit values, with 
only the signal yield and background parameters left free. 
Then, in each E*of^jj^ bin, the Tr'^ir/) correction factor is 



the ratio of the on-resonance minus luminosity-weighted 
off-resonance TT'~'{ri) yield to the luminosity-weighted MC 
Tr'^irj) yield. Systematic uncertainties from the fit are 
found by individually varying the fixed parameters in the 
on-resonance data fits, and also allowing for the MC fit- 
validation checks. The resulting correction factors and 
their uncertainties are shown in Tables [iTl and IIIII 

Correction factors to the BB MC predictions in 
100 MeV bins of E*, along with their uncertainties and 
correlations, are determined by applying the above fac- 
tors event-by-event to MC events passing the B — > Xg'f 
selection criteria. 



TABLE II: The tt" correction factors from ratios of data-to- 
MC fitted yields. The first and second sets of uncertainties 
are statistical and systematic, respectively. 



7r° CM Energy (GeV) 



Correction Factor 



1.4 to 1.6 
1.6 to 1.8 
1.8 to 2.0 
2.0 to 2.2 
2.2 to 2.4 
2.4 to 3.0 



0.959 ±0.006 ±0.013 
0.933 ± 0.009 ± 0.012 
0.990 ± 0.012 ± 0.031 
0.992 ±0.016 ±0.013 
0.899 ± 0.035 ± 0.018 
1.489 ± 0.259 ± 0.076 
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FIG. 6: Fits to spectra of 77 combinations per interval of 77 mass for tt" on-resonance data, in bins of E*o: (a) 1.4 to 1.6 GeV, 
(b) 1.6 to 1.8 GeV, (c) 1.8 to 2.0 GeV, (d) 2.0 to 2.2 GeV, (e) 2.2 to 2.4 GeV, and (f) 2.4 to 3.0 GeV. For each bin, the top plot 
shows the data (points), the total fit (upper curve) and the signal component of the fit (lower curve); the bottom plot shows 
the residuals, defined as (data - fit)/(data uncertainty). For the first few bins the signal shape doesn't precisely reproduce the 
center of the peak, an effect also seen in fits to MC 7r° signal only, but this does not affect the integrated signal yield. 



2. Additional corrections for low-energy photon efficiency 

While the procedure described above accounts for 
data-MC differences in the produced ir^ and rj yields af- 
ter the full selection, including the efficiencies for lepton 
tagging and for detecting the high-energy photon, it does 



TABLE III: The 77 correction factors from ratios of data-to- 
MC fitted yields. The first and second sets of uncertainties 
are statistical and systematic, respectively. 



■q CM Energy (GeV) 


Correction Factor 


1.5 to 1.7 


0.948 ± 0.029 ± 0.034 


1.7 to 1.9 


0.744 ± 0.026 ± 0.029 


1.9 to 2.2 


0.654 ±0.024 ±0.017 


2.2 to 2.6 


0.864 ± 0.049 ± 0.027 



not properly account for data-MC differences in the de- 
tection efficiency for the lovif-energy photon from a i:^ or 
77 decay. This is because the fits to the samples studied 
above count events in which that photon is detected and 
forms a 77 pair in the T^^irf) mass peak, whereas in the 
B — >■ Xs7 analysis a. B ^ XTi^{rf) background event is 
accepted if the low-energy photon is not found (or forms 
a reconstructed 77 pair mass outside the veto window). 
Thus the procedure corrects in the wrong direction for 
data-MC differences in low-energy photon detection effi- 
ciency. 

Correcting for low-energy photon efficiency is another 
multi-step process. First, BaBAR measurements of 7r° de- 
tection efficiency are taken from studies of the initial- 
state radiation (ISR) process e"^e^ — > u)^ with lj — > 
7r+7r~7r°. Here the precise knowledge of the beam en- 
ergies and the measured charged pions and high-energy 
ISR photon allow the four- momentum of the tt" to be pre- 
dicted. The measured efficiency difference between data 
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FIG. 7: Fits to spectra of 77 combinations per interval of 77 mass for 77 on-resonance data, in bins of E*: (a) 1.5 to 1.7 GeV, 
(b) 1.7 to 1.9 GeV, (c) 1.9 to 2.2 GeV, and (d) 2.2 to 2.6 GeV. For each bin, the top plot shows the data (points), the total 
fit (upper curve) and the signal component of the fit (lower curve); the bottom plot shows the residuals, defined as (data - 
fit)/(data uncertainty). 



and MC events is adjusted to match the tt" CM-frame 
momentum distributions for B — > Xtt" background in 
the B — >■ analysis and for the inclusive- tt" stud- 

ies described in section IVI B II The result is a data- 
MC fractional efficiency difference of (—4.1 ± 0.7)% for 
the B — > Xs^ selection and (—3.5 ± 0.6)% in the in- 
clusive 7r° studies. Part of these data-MC efficiency dif- 
ferences are accounted for by a data-MC difference of 
(— 1.15±0.65)% measured for the high-energy photon, as 
detailed in Sec. IX A II below. Subtracting this and com- 
bining the errors in quadrature leaves (—2.95 ± 0.95)% 
{B Xsi selection) and (-2.35 ± 0.9)% {B Xtt° 
selection) as due to the low-energy photon. 

Finally, the B Xn'^ samples are used to separately 
study the roughly 25% of low-energy photons in the cur- 
rent measurements that have laboratory-frame energies 
below 80MeV. This is necessary because, in order to 
suppress backgrounds, the ISR analysis effectively cov- 
ers cosines of the 7r° helicity angle (which equals the 
decay energy asymmetry) only up to about 0.9. Be- 
cause of this, low-energy photons below « 80MeV are 



not adequately represented in the ISR analysis. While 
the B — > Xtt^ samples cannot provide the overall low- 
energy photon correction, they can be used to determine 
the relative correction between photons below 80MeV 
and those above BOMeV, by separately applying the tt" 
mass-spectrum fitting technique for tt" mass combina- 
tions involving low-energy photons in these two regions. 
An additional data-MC fractional efficiency correction of 
(— 3.6±l.l)%is derived for only those tt" decays involving 
these photons below 80 MeV. There is no corresponding 
effect for 77's, where the minimum photon energy is al- 
ways at least 75 MeV. 

To determine the effect of these low-energy-photon ef- 
ficiency differences on the analysis, the BB simulation is 
rerun with the specified fractions of low-energy photons 
from 7r'^(?7) decays discarded. The result is an additional 
factor of 1.105±0.029 for the 7r° component of BB back- 
ground in the B ^ Xgj analysis, and 1.041 ± 0.015 for 
the ry component. The 7r° and ij errors are mostly cor- 
related. These factors multiply those obtained from the 
inclusive ir'^i'rj) data/MC yield comparisons. 
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C. Other Meson Decays 

Radiative decays of inclusively produced uj (in the n'^j 
mode) and 77' (in various decay modes) can lead to high- 
energy photons not already accounted for among the in- 
clusive 7r°'s. As seen in Table HI these contribute sev- 
eral percent of the simulated BB background. We have 
studied inclusive uj and rj' production in T(4S') events. 
Correction factors are determined in bins of CM-frame 
meson momentum (p* or p*,) as the ratios of measured 
inclusive branching fractions to the values used in the 
MC simulation for the current analysis. 

The Lo measurements cover the range from to 
2.25GeV/c in 0.25GeV/c bins. Correction factors range 
from 0.7 to 1.3, with uncertainties averaging 0.17. 

Results for rj' are divided into two regions of reduced 

momentum, Xn> = P'^' / \J E^eam^ ~ For 0.39 < Xr,' < 
0.52, direct B — >■ Xrj' decays are dominant, and the cor- 
rection factor is 1.86±0.61. For 0.10 < Xjj' < 0.39, decays 
via an intermediate charm-meson state are dominant, 
and the correction factor is 0.35 ±0.19. The first range is 
most important in the signal region for B — Xs'y, while 
the second range is most important in the BB control 
region. 

Both w and rj' corrections are applied event-by-event in 
the BB simulation in order to obtain correction factors 
in EZ bins. 



D. Electron Backgrounds 

Electrons and positrons contribute to the photon back- 
ground in two ways (see Table U). First, there are 
events in which the reconstructed photon is from hard 
bremsstrahlung from an interacting with the material 
in the inner portion of the BaBAR detector (beam pipe, 
SVT, and material between the SVT and the active area 
of the DCH). Second, there are events in which the recon- 
structed photon is faked by an electron due to a failure to 
reconstruct a track or to match a track to the calorimeter 
energy deposit. The primary source of the in both of 
these categories is semileptonic B decay. 

The bremsstrahlung process is reliably simulated by 
GEANT4, so there is no correction to the simulation for 
this background. But a 3% systematic error is assigned 
based on the precision with which the amount of detector 
material has been measured. 

The misreconstructed electron background is measured 
using a tag and probe method with B — >■ XJ/ip{J/^ — 
e^e^) data. This sample closely models the particle mul- 
tiplicity in i? — >■ Xgj events. The J/ip in this decay mode 
is normally reconstructed by requiring two electrons with 
tracks associated with EMC clusters. If the track is mis- 
reconstructed there will still be a cluster but without an 
associated track. In this case the J/tp is reconstructed 
from this unassociated cluster along with the other elec- 
tron, which has a track matched to a cluster. Because 



either of the two leptons could have a misreconstructed 
track, the track inefficiency may be measured as 

7V(J/i/'(eClus,eTrk)) 
^ " ~ 27V(J/V'(eTrk,eTrk)) + Ar( J/V'(eClus, eTrk)) ' 

(7) 

where iV( J/?/>(eClus, eTrk)) and iV( J/V'(eTrk, eTrk)) are 
the numbers of J/ip — >■ e"'"e~ events with one and two 
reconstructed tracks, respectively. These yields are ex- 
tracted from fits to distributions of e'^ e~ invariant mass, 
computed from the four-momenta of the track found for 
one lepton (the "tag") and the EMC cluster for the other. 
The value of 1 — e is compared between data and MC 
samples to derive a correction factor for the simulation. 
There is a large combinatoric background in the one- 
track (eClus, eTrk) sample due to actual photons. How- 
ever, when an electron track has been misreconstructed 
there are still a number of DCH hits around the tra- 
jectory from the vertex to the EMC cluster. The back- 
ground is significantly reduced by requiring a minimum 
number of 20 hits in a road of 1-cm radius around this 
trajectory. 

Figure [5] shows an example of fits to the e~^e~ mass 
combinations corresponding to the numerator and de- 
nominator of Eq. ([7]) for data and MC simulation. The 
mass is computed from the track associated with a tag 
lepton and the EMC cluster associated with the other lep- 
ton; hence (eTrk, eTrk) combinations are entered twice, 
with different masses, once for each tag. The simulation 
underestimates the fraction of misreconstructed tracks 
by a factor of 1.57± 0.27(stat) ± 0.22(syst), where "stat" 
and "syst" denote the statistical and systematic uncer- 
tainties, respectively. The systematic error comes pre- 
dominantly from uncertainties in the line shape assumed 
in the invariant mass fit and from varying the require- 
ments on the road width and the number of DCH hits. 
Consequently, the MC estimate of the BB background 
to i? — > Xs"f from misreconstructed electrons is increased 
by a factor of 1.57 ± 0.35. 



E. Antineutrons 

The only significant hadron background to high-energy 
photons is from antineutrons, which have a neutral signa- 
ture and can, by annihilating in or just before the EMC, 
deposit a large amount of energy. A large fraction of such 
background is removed by the requirement on maximum 
lateral moment fSec. lIVB]) . There are two sources of po- 
tential bias in the predicted yield: the inclusive B — )■ Xn 
branching fraction and n momentum spectrum in the 
event simulation, and the GEANT4 simulation of the 
deposited energy and its distribution in the EMC. Be- 
cause it is not possible to identify or measure the four- 
momentum of an n in the BaBAR detector, there are no 
control samples of n's available to study these effects. 
Hence estimates of their size have been based on com- 
parison of data to simulated events involving B decays 
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FIG. 8: Invariant e^e~ mass distributions for J/i/) — >■ e'''e~ electron-inefficiency studies. Plots (a) and (b) are for MC samples; 
plots (c) and (d) are for data samples. Plots (a) and (c) includes all (eClus, eTrk) and (eTrk, eTrk) pairs, while plots (b) and 
(d) include only the (eClus, eTrk) pairs. Dotted, dashed and solid curves show, respectively, the J/^ signal and background 
components of each fit, and their sum. The numbers of fitted J/^p signal events are 156327 ± 550, 313 ± 27, 35825 ± 281 and 
109 ± 18 for plots (a), (b), (c) and (d), respectively. 



to p's. 

The inclusive n production spectrum in the BB sim- 
ulation is corrected by the ratio of a measured inclu- 
sive p spectrum to its corresponding simulation. Correc- 
tion factors are applied as a function of CM-frame an- 
tibaryon momentum. They are close to 1.0 at momenta 
above about 0.9GeV/c, increasing to 2.0 for momenta 
from 0.5 GeV/c down to the lowest measured momentum 
of about 0.3GeV/c. Uncertainties are typically 8% to 
12%. Below 0.3GeV/c, a factor of 2.0 is assigned, with 
a larger uncertainty. In addition, while the production 
of p and n from direct B decays is related via isospin 
conservation, many of the antibaryons arise from decays 
of A's or hyperons, which would require separate correc- 
tion factors. An additional uncertainty of 3% accounts 
for differences in fractions of direct B vs. A vs. hyperon 
parentage of n and p. 

Control samples of p's from the decay of A's are used to 
compare data and MC EMC response to p's as a function 
of laboratory-frame p momentum. Most p's are rejected 
by imposing the same upper limit on the lateral moment 
of their EMC energy deposition pattern as used in the 
B — >■ Xg'y photon selection. Correction factors are de- 
termined in bins of laboratory-frame momentum pp vs. 



xemc, where 



^T^ + ml +mp 

is that fraction of the total energy from annihilation on 
a nucleon that is deposited in the EMC. Corrections are 
computed as the ratio of data to MC probability distri- 
bution functions (PDFs) for xemc > 0-5, the only re- 
gion that can yield an apparent E* above 1.53 GeV. The 
primary data-MC differences result from the larger av- 
erage lateral moment in data, an effect accentuated by 
restricting the lateral moment to low (photon-like) val- 
ues. Hence the data have a considerably smaller propor- 
tion of p's satisfying the selection than is the case for 
MC events. This inaccuracy of the simulation, and to 
a lesser extent an overestimate of the energy deposit it- 
self, increases with increasing XemCi hence the data/MC 
correction factor becomes small as xemc increases. 

Figure [5] illustrates the p correction factors in this two- 
dimensional space. Note that in inclusive B decays there 
are relatively few antibaryons with laboratory-frame mo- 
mentum above 1.5GeV/c. 

There are several difficulties in applying these results 
to antineutrons. First, because of energy loss in the beam 
pipe and inner detector components, p's do not provide 
useful energy-deposit information for laboratory-frame 
momenta below 0.5 GeV/c, whereas one-third of Ti's from 
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FIG. 9: Data/MC correction factors for |5's in three illustra- 
tive intervals of laboratory momentum, vs. a;EMC (see Eq. [8|, 
from A control samples, including the same maximum lateral 
moment criterion used for B — )■ Xs'y analysis. 



B decay have momenta below 0.5 GeV/c. A constant ex- 
trapolation of correction factors to lower momenta is as- 
sumed, with a systematic uncertainty set by including 
an additional factor of 1/2. Second, p's enter the EMC 
crystals at a larger angle of incidence than do n's, be- 
cause of the magnetic field, resulting in a larger lateral 
moment for p's with px < 0.7GeV/c than for n's. A sys- 
tematic uncertainty is assigned by increasing correction 
factors with laboratory pp in this region to their values 
at just- higher, unaffected, pp. Third, because of a mis- 
take in the version of GEANT4 implemented in BaBAR, 
simulated p's which stop before annihilating do not then 
annihilate. This has been dealt with by increasing the 
MC PDFs (decreasing the correction factors) according 
to the fraction of p's which annihilate for a given momen- 
tum. Half of this correction is adopted as its systematic 
uncertainty. 

Correction factors to the simulated n background in 
E* bins are computed by applying event-by-event cor- 
rections for both the branching fraction and the EMC 
response. Systematic uncertainties are obtained by re- 
doing this for each of the systematic changes outlined 
above. The resulting correction factors vary from about 
0.4 to 0.04 as E* increases from 1.53 to 2.8 GeV, with 
uncertainties ranging from 1/4 to 1/2 of the correction 
factors. 



F. Final-State Radiation 

Final-state radiation, most importantly from leptons, 
is incorporated into the BB background simulation with 
PHOTOS 53]. The contribution is labeled as having B 
parentage in Table HI No correction is applied for this 
small component. Radiation from light quarks during 



the hadronization process is not incorporated into the 
simulation. However, this contribution was computed 
for the previous B ^ Xgj analysis [i^l, where a photon 
spectrum based on the calculation in reference [54] was 
passed through the detector simulation and selection cri- 
teria. This contribution was found to be less than 0.3%. 



G. Semileptonic Branching Fraction 

The dominant source of tagging leptons above the min- 
imum required momenta (Sec. IIV Cp in both signal and 
BB background events, and also of electrons that fake 
high-energy photons (Table H]) is the semileptonic decay 
of B mesons. The MC simulation models B semilep- 
tonic decays as a sum of exclusive processes. But this 
sum does not accurately reproduce inclusive measure- 
ments of semileptonic decays [sst - A BaBAR inclusive 
electron measurement [sgI [57| is used to renormalize the 
simulated branching fractions as a function of CM-frame 
lepton momentum Figure [TU] shows the data and MC 
points and their ratio. Correction factors are applied 
based on the polynomial fit. For most leptons relevant 
to this analysis the correction is larger than unity. 

This correction enters in two places in the analysis. 
First, it affects tagging efficiency. By integrating over all 
lepton tags in events passing selection criteria, a correc- 
tion factor of 1.047 ± 0.013 is obtained, for B -> Xg'f 
signal events, while for the BB MC sample the factor is 
1.051±0.013. This correction is independent of i?* . How- 
ever, the procedure for normalizing the tt" and rj back- 
ground components to data implicitly takes this into ac- 
count. Hence the correction is applied only to other BB 
components. (However, the corrections given in Table U 
for these components are derived before applying this ad- 
ditional semileptonic correction.) 

In addition to its effect on lepton tagging, the semilep- 
tonic correction affects the two backgrounds in which an 
fakes a high-energy photon. The corrections (which 
are included in Table H] along with the corrections 
described in Sec. IVID|) depend upon E*; their average 
value for both backgrounds is 1.058 ± 0.013. 

The two effects are taken to be fully correlated in com- 
puting their contribution to the overall BB yield uncer- 
tainty. 

H. Overall BB corrections 

The above subsections describe corrections by BB 
component for all but a few percent of the predicted 
makeup of the BB background, as summarized in Ta- 
ble m Several other small categories (e.g., "J/ip" and 
"Other" ) are left as-predicted. Finally, several small cor- 
rections computed in the context of signal efficiency are 
also applicable to BB backgrounds: a high-energy pho- 
ton efficiency correction of 0.9885 ± 0.0065 (Sec. IX A ip . 
and a correction of 0.989 ± 0.004 for lepton identification 
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FIG. 10: _B-meson semileptonic partial branching fraction vs. 
CM-frame lepton moment um, averaged over B charge states. 
Top: BABAR measurement f57| (filled black circles) and values 
in the BB MC simulation (open circles). Bottom: data/MC 
ratios, and results of a second-order polynomial fit from 1.0 
to 2.2 GeV. The dashed curves show the la error band. 



efficiency in a multi-particle environment (Sec. IX A3|) . 
Like the semileptonic tag correction, these need only be 
applied to the 20 to 25% of B backgrounds other than 
TT*^ or 7y, and hence are small effects. More significant is 
a global factor of 0.991 ± 0.004 from different probabil- 
ities between MC and data events of the 7r° veto being 
activated by a background photon. Uncertainties also 
include a small contribution from BB MC statistics. 

The BB corrections described above are applied to 
each component and for each 100 MeV bin of E*. Cor- 
relations between bins, due both to £'*-independent cor- 
rections and to corrections dependent upon parent en- 
ergies, are tracked, resulting in a table of corrected BB 
yields and a correlation matrix. This information is used 
to compute the results presented in the next section. 
The largest systematic uncertainties on the BB yields 
are those due to the low-energy photon correction to the 
TT^ and r] components, with uncertainties in the no-track 
electron component and the tt'^ inclusive spectrum next 
most significant. 
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VII. SIGNAL YIELDS AND VALIDATION OF 
BACKGROUND ESTIMATION 

Figure [11] shows the photon energy spectrum in data 
after background subtraction. Table ITVl gives the signal 
yields and background estimations in bins of E*. The 
bin-to-bin correlations between the errors on the signal 
yields are given in Table|Vl The continuum background is 
estimated with off-resonance data, while the BB back- 
ground is estimated from MC simulation, with all the 
corrections described in Sec. |VT] applied. 

To validate the background estimation, two control 
regions are set aside in the photon spectrum. In the 
upper control region (2.9 < E* < 3.5 GeV), the event 
yield after subtracting continuum and BB backgrounds 
is —100 ± 138 (stat) ± 14 (syst) events, where the statis- 
tical uncertainty results from off-resonance subtraction. 
The systematic error is from the uncertainty of out-of- 
time Bhabha-scattering events in BB background (see 
Table H] caption). This is consistent with the expectation 
of zero events. 

In the lower control region (1.53 < E* < 1.8 GeV), 
there remain 1174 ± 272 (stat) ± 828 (syst) events after 
background subtraction. The errors in the BB estimates 
in these E* bins are highly correlated; these correlations 
have been included when computing the control-region 
systematic error. The agreement with zero in this region 
is at the 1.4a level, assuming no signal events. However 
this energy region contains a few hundred signal events, 
with the exact number depending on the assumed signal 
model. For example, using predictions based on the ki- 
netic and shape function schemes with parameters close 
to HFAG's world-average values [11], on average about 
275 signal events would be expected in the lower control 
region. Allowing for this, the data-background difference 
is reduced to the I.Oct level. 
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TABLE IV; The event yields in bins of E*. The continuum background is estimated from off-resonance data normalized to 
on-resonance luminosity. The BB background is estimated using Monte Carlo simulation, corrected as described in Sec. IVII 
The extracted signal yield is computed by subtracting the continuum and BB backgrounds from the on-resonance data yield. 
It is quoted with statistical uncertainties (from on-resonance minus off-resonance subtraction) and BB systematics. The first 
three rows show updated numbers for the control region. The last set of rows show yields in wide E* bins, taking into account 
the correlations of BB backgrounds between 100-MeV bins. 



E* (GeV) 


On-resonance 
Data 


Continuum 
Background 


BB 
Background 


Signal 
Yield 


1.53 to 1.60 


11869 ± 109 


1319 ± 112 


10232 ± 275 


318 ± 156 ± 275 


1.60 to 1.70 


13531 ± 116 


1327 ± 113 


11497 ±316 


706 ± 162 ± 316 


1.70 to 1.80 


10366 ± 102 


1371 ± 115 


8846 ± 252 


150 ± 153 ± 252 


1.80 to 1.90 


8054 ± 90 


1118 ± 105 


6511 ± 195 


426 ± 138 ± 195 


1.90 to 2.00 


6083 ± 78 


885 ± 93 


4732 ± 139 


466 ± 121 ± 139 


2.00 to 2.10 


4429 ± 67 


717 ±82 


3165 ± 91 


548 ± 106 ± 91 


2.10 to 2.20 


3124 ± 56 


659 ± 80 


1743 ± 56 


722 ±98 ±56 


2.20 to 2.30 


2465 ± 50 


603 ± 77 


757 ± 33 


1105 ±91 ±33 


2.30 to 2.40 


1977 ± 45 


639 ± 79 


314 ± 20 


1024 ±90 ±20 


2.40 to 2.50 


1712 ± 41 


537 ± 73 


152 ± 19 


1024 ±84 ±19 


2.50 to 2.60 


1225 ± 35 


499 ± 71 


67 ±9 


659 ±79 ±9 


2.60 to 2.70 


795 ± 28 


328 ± 55 


32 ±7 


435 ±62 ±7 


2.70 to 2.80 


457 ± 21 


404 ± 62 


18 ±3 


35 ±66 ±3 


2.80 to 2.90 


410 ± 20 


310 ±55 


9±4 


91 ± 59 ± 4 


2.90 to 3.00 


370 ± 19 


292 ± 52 


8±4 


71 ± 55 ± 4 


3.00 to 3.10 


298 ± 17 


335 ± 56 


6±3 


-44 ± 59 ± 3 


3.10 to 3.20 


305 ± 18 


396 ± 61 


5±3 


-96 ± 64 ± 3 


3.20 to 3.30 


279 ± 17 


273 ± 51 


6±2 


± 54 ± 2 


3.30 to 3.40 


252 ± 16 


318 ± 56 


3±2 


-69 ± 58 ± 1 


3.40 to 3.50 


222 ± 15 


182 ± 42 


3± 1 


38 ± 44 ± 1 


1.80 to 2.80 


30321 ± 174 


6387 ± 249 


17490 ± 496 


6444 ± 304 ± 496 


1.90 to 2.80 


22267 ± 149 


5270 ± 226 


10980 ± 313 


6018 ±271 ±313 


2.00 to 2.80 


16184 ± 127 


4385 ± 206 


6248 ± 187 


5551 ± 242 ± 187 


2.10 to 2.80 


11755 ± 108 


3669 ± 189 


3083 ± 110 


5004 ±218 ±110 



VIII. OBTAINING PHYSICS RESULTS: AN 
OVERVIEW 

Three physics results are extracted from the measured 
signal yield: 

• the CP asymmetry, Acp {B Xs+dl), 

• the inclusive branching fraction, B{B — > Xs^) (for 
several wide ranges of true E,y in the i?-meson rest 
frame), and 

• the true spectral shape and energy moments for 
B — > (in both the CM frame and the B frame). 

The presence of new physics beyond the SM can affect 
the branching fraction and Acp- The spectral shape, 
however, depends only on the dynamics of the b quark 
within the B meson; it is independent of any new physics 
contributions. Three different approaches are optimal for 
the three physics results. 



The branching fraction and spectral shape measure- 
ments require corrections for efficiency. The partial 
branching fraction for signal in any range of measured 
photon energy E* is obtained from the signal yield S in 
that same range by 

B{B ^ X,+al) ^ ,r^^ , (9) 

where Sgig is the signal efficiency for that range and N^-g 

is the number of BB events in the on-resonance data set 
before event selection. B{B Xgj) is obtained from this 
by removing the small constant fraction contributed by 
B — > X(i7. Applying Eq. ^ brings in additional system- 
atic uncertainties related to the efficiency and to N^-g. 
The inclusive branching fraction and spectral shape mea- 
surements are made in terms of reconstructed E* in the 
CM frame, while theoretical predictions are made for true 
photon energy E^ in the B frame. These differ due to res- 
olution and Doppler smearing. The measurements must 
be converted to corresponding measurements in terms of 
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TABLE V: The correlation matrix for the signal yield errors from Table llVl in 100-MeV bins of E* . Systematic {BB background) 
and statistical contributions are included. Rows and columns are labeled by the value of E'^ at the lower edge of the bin. 



e; (GeV) 


1.53 


1.6 


1.7 


1.8 


1.9 


2.0 


2.1 


2.2 


2.3 


2.4 


2.5 


2.6 


2.7 


1.53 


1.00 


0.75 


0.71 


0.65 


0.58 


0.46 


0.32 


0.16 


0.07 


0.03 


0.01 


0.01 


0.00 


1.6 




1.00 


0.74 


0.68 


0.61 


0.48 


0.33 


0.17 


0.07 


0.03 


0.02 


0.01 


0.00 


1.7 






1.00 


0.67 


0.60 


0.47 


0.33 


0.17 


0.08 


0.03 


0.02 


0.01 


0.00 


1.8 








1.00 


0.58 


0.46 


0.32 


0.17 


0.08 


0.04 


0.02 


0.01 


0.00 
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0.01 
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0.00 


0.00 
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0.00 
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true E* or , in order to allow for detector-independent 
comparisons. 




E/ (GeV) 

FIG. 12: Comparison of data spectrum in reconstructed E* 
to the predictions of three models, each normalized for best 
agreement with the data above 1.8 GeV, based on (in- 
cluding bin-to-bin correlations). The solid histogram is for 
a shape-function-scheme model with mj, — 4.51 GeV and 
fj,^ = 0.46 GeV^ which resembles the data in this range. The 
dot-dash histogram is for a kinetic scheme model with param- 
eters mb = 4.60 GeV and fil = 0.45 GeV^, close to the HFAG 
world average. The dotted histogram is for a shape-function- 
scheme model with rub = 4.40 GeV and fil = 0.52 GeV^. The 
minimum values of are 6.7, 13.4 and 19.6, respectively. 



Efficiency factors and also the transformation from one 
definition of photon energy to another depend upon the 
choice of signal model, i.e., on the values of the HQET 
parameters in the kinetic or shape function scheme. 
HFAG [1^ has extracted world-average values of these 
parameters by combined fits to measurements of i? ^ 
Xc^v decays and previous B — > Xs^ measurements. For 



the present inclusive branching fraction measurement the 
range of models considered is based on these HFAG cen- 
tral values and errors. On the other hand, for the spec- 
trum measurement such a restriction would prejudice the 
results; the range of models considered must instead be 
driven by the data. Put another way, for the branching 
fraction measurement the MC model plays a subsidiary 
role, used only to estimate the efficiency and transforma- 
tion factors, so it makes sense to use the best available 
information to constrain the model, while for the spec- 
trum the model is itself the object of the measurement. 
With these procedures, the model-dependence uncertain- 
ties for both measurements are small compared to the 
combined statistical and systematic uncertainties. 

In Fig. [T^l predictions of three models are superim- 
poses on the measured data. The first resembles the data 
for measured E* above 1.8 GeV. The second, which has 
HQET parameters very close to the HFAG world-average 
values in the kinetic scheme, is about one standard devi- 
ation ("Icr") below the data in the first few energy bins 
above 1.8 GeV, where BB background is large. The third 
is somewhat more than Icr above the data in this re- 
gion. Differences between data and a particular model 
may be due either to the model being an incorrect de- 
scription or to systematic fluctuations in the BB back- 
ground contribution. This recognition is a key element of 
the approaches used to measure both the branching frac- 
tion (see below) and the shape of the true energy spectra 
(Sec.KEI. 

Branching fraction results are determined for E* > 
1.8 GeV. The branching fraction is computed by apply- 
ing Eq. dH) to a single wide bin, e.g., 1.8 < E* < 2.8 GeV, 
using the average efficiency Cgig computed for an HFAG- 
based model. If l/csig factors were instead applied in 
100-MeV bins, the smaller values of Cgig at low energies 
(Fig. [5]) would amplify the larger systematic uncertainties 
on the event yield in this region as well as any data-model 
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differences, in effect translating possible background fluc- 
tuations into a larger branching fraction bias. Because 
of the energy-dependent Csig, statistical precision also im- 
proves with fewer bins. Note also that the model depen- 
dence of the branching fraction computed using 100-MeV 
bins is comparable to that for a wide bin. Thus, overall, 
the wide-bin approach is both more accurate and more 
precise. Full details are in Sec. [Xj 

In contrast, the spectral shape must be determined by 
applying Eq. IHlin each 100-MeV bin of reconstructed E*. 
This is the first step of a four-step unfolding procedure, 
detailed in Sec. |Xll leading to the true photon energy 
spectra in the CM and i?-meson rest frames. Each model 
shown in Fig. [T^] is used in all four steps, to obtain the 
measured spectrum and its model dependence. Energy 
moments and their correlations are computed from the 
unfolded spectra. This information is a needed input 
to the HFAG fitting procedure, and may facilitate other 
potential comparisons with theory. 

The effects of efficiency and smearing cancel in the ex- 
traction of Acp. A raw asymmetry is thus directly com- 
puted from the measured yields vs. E* , using the lepton 
charge to tag B vs. B mesons. Systematic corrections 
and uncertainties arise only from possible charge depen- 
dence of the efficiencies (which would be a bias), as well 
as from mistagging (which dilutes the asymmetry). The 
full Acp analysis procedure is described in Sec. |lXl 



IX. MEASUREMENT OF DIRECT CP 
ASYMMETRY 

The direct CP asymmetry, Acp {B Xg+dj) is mea- 
sured by dividing the signal sample into B and B decays 
according to the charge of the lepton tag and computing 

A^^iB ^ X,+,7) = ^ll^- . (10) 

where iV+'^^^ are the positively (negatively) tagged signal 
yields. The asymmetry must be corrected for the dilution 
due to the mistag fraction w: 

AcpiB ^ X.,+dl) = Y^Agr (5 ^ X,+rf7) . (11) 

As can be seen in Fig. [2jb) and Table IIVI the BB back- 
ground decreases at higher photon energies. It was deter- 
mined (prior to looking at the data) that restricting the 
Acp signal region to 2.1 < E* < 2.8 GeV optimizes the 
statistical precision, and also the total precision includ- 
ing the uncertainty on the BB background asymmetry 
described below. Other systematic uncertainties on Acp 
have negligible variation with E*. The theoretical SM 
prediction of a near-zero asymmetry is not affected for a 
minimum-energy requirement of 2.1 GeV (ol. [ssj. All of 
the other selection requirements (Sec. IIVP were found to 
be optimal also for the Acp measurement. 



The tagged signal yields are N+ = 2620 ± 158(stat) 
and N~ — 2389 ± 151(stat), giving an asymmetry of 

Ag^''"(B -> Xs+dj) = 0.046 ± 0.044 . (12) 

To correct for dilution we compute the mistag fraction 

1^ = ^ + ti^cascadc + WmisID • (13) 

The largest contribution is from B"-B° oscillation, with 
mixing probability Xd = 0.1863 ± 0.0024 [55]; the factor 
of 1/2 accounts for the B^ mesons, which do not oscil- 
late. Smaller contributions are Wcascado = 0.0328±0.0035, 
the fraction of events with wrong-sign leptons from the 
B decay chain, and uij^isiD = 0.0073 ± 0.0037, the 
mistag fraction due to misidentification of hadrons as lep- 
tons (almost entirely in muon rather than electron tags). 
Their values are taken from the MC simulation averag- 
ing over electron and muon tags. An additional uncer- 
tainty in uj arises because our MC simulation assumes 
B{T{4:S) B"W) = 0.50 which leads to the factor of 
1/2 in the first term of Eq. (IT^ . The measured value 
is S(r(4S') ^ B^B°) = 0.484 ± 0.006 [H so we take 
as a systematic the difference between the measured and 
assumed values, Acj = 0.016xd- This and the errors on 
Xd, ^cascade and WmisiD are added in quadrature to give 
a; = 0.131 ±0.006. _ 

The uncertainty in the BB background estimation de- 
scribed in Sec. I VII cancels in the numerator of Eq. (fTO)) 
but not in the denominator, leading to an uncertainty in 
Aqp^ of 0.022. This uncertainty is combined with the 
uncertainty in uj to give a multiplicative systematic un- 
certainty on AcpiB X,+dl) of 0.029Acp. Table EJ 
summarizes all of the contributions to lj and to this un- 
certainty. 

The measured asymmetry could be biased if there were 
(a) an asymmetry in the BB background not modeled in 
the simulation or (b) a charge asymmetry in the lepton 
tag efficiency. To assess the potential bias due to BB 
subtraction the asymmetry of the BB data in the con- 
trol region 1.53 < E* < 1.8 GeV is measured in data to 
be A^^''^{B control) = 0.006 ± 0.009(stat). Interpreting 
this as a bias, it translates to a correction for the Acp 
signal region of AA^^^^B X,+dj) = -0.004 ± 0.006. 
The B — T' TT^X background sample described in Sec. lVI B] 
is used to confirm that there is no E* dependence to this 
correction in the signal region. Lepton charge tag asym- 
metries have been measured in e^e~ — >■ e"'"e~7, e~^e~ — >■ 
/i^7 and B K^*^J/i}){l+i^) events. No significant 
asymmetries are observed to a precision of 0.011, which 
is assigned as a systematic error on A^'i^{B — Xg+dl)- 
Table IVIII summarizes these additive systematic effects, 
showing a combined error in quadrature of 0.013. 

Since the SM prediction of Acp ~ depends upon 
cancellation of i? — > ^^7 and B — > Xd') asymmetries, a 
difference in their selection efficiencies could also cause a 
bias. We have used MC simulations with the same un- 
derlying model (KN with mb = 4.65 GeV/c^) to compare 
selection efficiencies following sq vs. dq hadronization. 
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TABLE VI: Contributions to Acp multiplicative systematic 
correction and error. 



Source 




A Acp /Acp 


B°B° oscillation 


(0.1863 ±0.0023)/2 




Fake lepton ID 


0.0073 ± 0.0037 




Cascade decays of B's 


0.0328 ± 0.0035 




B°W : B+B- = 1:1 


0.0000 ± 0.0030 




Total Lj 


0.133 ±0.0064 


0.018 


BB yield 




0.022 


Total uncertainty 




0.029 



Section IX Al describes systematic corrections and un- 
certainties affecting the efficiency Csig in Eq. and com- 
putes the total fractional systematic uncertainty on the 
branching fraction. The choice of the central values for 
Esig and a depend upon the choice of the signal model 
used in MC simulation. Section |XB] addresses this choice 
and determines the model-dependence uncertainty of the 
branching fraction. Section IX CI presents branching frac- 
tion results first in terms of measured E* , then presents 
the conversion to the branching fraction in the B-meson 
rest frame, along with associated uncertainties. 

A. Systematic Corrections and Uncertainties 



TABLE VII: Acp^^ additive systematic corrections and errors 



Source 


Correction (10"^) 


BB background 


-0.4 ±0.6 


Detection asymmetry 


0.0 ± 1.1 


Total 


-0.4 ± 1.3 



For E* > 2.1 GeV, the B -J> Xd^y efficiency is larger by 
a factor of 1.028 ± 0.014, so we conservatively assign a 
4.2% uncertainty (MC central value plus one standard 
deviation) in the yield oi B X^j events. Given the 
SM-predicted yields and asymmetries [l2|, that would 
change A'^p'^ by less than 0.0002, which is negligible. 

Finally the A^p^^{B — )• Xs+dl) is corrected for mistags 
and bias to give 

. (0.046 ± 0.044(stat)) - (0.004 ± 0.013) 

0.734(1 ±0.029) 
= 0.057 ± 0.060(stat) ±0.018(syst) , 

where the two systematic errors have been combined in 
quadrature. The result is consistent with no asymmetry. 



X. MEASUREMENT OF B{B Xsj) 

As discussed in Sec. IVIIII B{B — > ^^+^7) is mea- 
sured by applying Eq. to a single wide bin in mea- 
sured E* . Results are computed for three choices of en- 
ergy range: 1.8 to 2.8 GeV, 1.9 to 2.8 GeV and 2.0 to 
2.8 GeV. Note that egig here means the overall signal effi- 
ciency, including both acceptance and event-selection, as 
discussed in Sec. IIVFI A small adjustment, by a fac- 
tor a which is close to 1.0, converts each result to a 
branching fraction in the same range of the true Ej in 
the S-meson rest frame. This corrects for the effects 
of EMC resolution and Doppler smearing. Finally, the 
factor 1/(1 ± (|Vtd|/|Vts|)^) is applied to account for the 
contribution of B — >■ X^'j events, yielding a branching 
fraction for B — )■ only. 



Each of the factors in Eq. (jH]) can contribute to the 
uncertainty in the branching fraction. The signal yield 
has contributions from statistics (of the on-peak and off- 
peak data yields) and from the systematics of the BB 
background subtraction. The number of produced T{AS) 
events, N^-^, has a systematic uncertainty of 1.1%. The 
focus here is on the systematic uncertainty of the re- 
maining factor, the signal efficiency Csig- For each event- 
selection criterion, an efficiency is computed using MC 
simulation. But the actual efficiency in data may dif- 
fer from that in the simulation. Systematic corrections 
are determined by comparing data to MC events for var- 
ious control samples; the precision of each comparison 
provides a systematic uncertainty. A summary of these 
corrections and uncertainties is presented in Table IVIIII 

1. Systematics of the High-Energy Photon Selection 

Two dedicated studies of high-energy photon detection 
efficiency have been done using /i/i7 ISR events. These 
events are overconstrained, so the measured ^+ and jj," 
tracks, along with known beam kinematics, can be used 
in a IC fit to predict the three-momentum of the photon. 
Naively one would look for a detected photon "close" 
to this predicted photon, and the data/MC correction 
would be the ratio of the probability of finding such a 
photon in fi/ij data events to the probability of finding 
one in /i/i7 MC events. However, this is complicated by 
the effects of EMC resolution and by the possibility that 
the likelihood of photon conversion (in detector material) 
isn't accurately simulated. The earlier and more recent 
of the studies took rather different approaches to 
these issues. The first applied acceptance criteria (par- 
ticularly a minimum energy) to detected photons, and 
folded EMC resolution into the predicted photon prop- 
erties before making the same cuts. That study also in 
effect measured the conversion fraction, separately for 
data and MC samples. The second study did not use 
acceptance cuts for the detected photon, instead loosely 
matching its parameters to those of the predicted pho- 
ton, and used an electron veto to suppress photon conver- 
sions. Results of the two studies are in good agreement. 
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with the data/MC efficiency corrections differing by 0.3% 
when weighted by the B — ^ Xsj photon polar angle dis- 
tribution. Systematic uncertainties on these corrections 
are 0.65% and 0.55%. The two correction factors are 
averaged, giving 0.9885, and an uncertainty of 0.65% is 
assigned. The assigned correction and uncertainty are 
independent of the photon energy E*. Hence they af- 
fect the branching fraction, but not the spectral shape or 
energy moments. 

The samples from data and MC simulation are 
also used to assess the photon energy scale and resolu- 
tion, by comparing the distributions for data and MC 
events of the ratio of detected to predicted photon en- 
ergy. For the energy scale, the energy balance in the 
decay K*{892){-^ TT+n~)j is also used. After 

small energy-scale adjustments already included in event 
reconstruction, both processes show no remaining bias 
for either MC or data events, with a conservative uncer- 
tainty of 0.3%. For photon energy resolution, inclusion of 
an additional 1% energy smearing of MC photons brings 
the fifij ratio distribution into good agreement with that 
for fj,fj,^ data. This is taken as a systematic uncertainty. 
The energy-scale and resolution effects translate into the 
small uncertainties on the inclusive branching fraction 
shown in Table |VllIl 

Lastly, the samples are used to assess shower 

shape, in particular the efficiency of the selection cut on 
lateral moment. After a small adjustment of the simu- 
lated lateral moment, there is good agreement between 
MC and data efficiencies of this selection, with the un- 
certainty given in Table IVIIII 

The high-energy photon efficiency is calibrated using 
the low-multiplicity events, but could also be af- 
fected by the hadronic-event environment in BB events 
(including signal). The requirement that the high-energy 
photon be isolated from any other EMC energy deposi- 
tion by at least 25 cm is meant to reduce data-MC effi- 
ciency differences. The systematic uncertainty of 2% is 
estimated by embedding high-energy photon signatures 
into hadronic events, separately for data and MC sam- 
ples, and determining the fractions of events passing the 
isolation requirement. 



2. Systemattcs of the vr and rj Vetoes 

The 7r° and rj vetoes can remove events not only if the 
high-energy photon originates from an actual tt" or ?/, 
but also if there is a random ( "background" ) photon with 
which the high-energy photon forms a 77 invariant mass 
combination lying inside one of the veto windows. The 
efficiencies of the vetoes for simulated events can differ 
from those for data if the number of background photons 
in simulation differs from data. Off-resonance-subtracted 
data and BB MC events are compared for high-energy 
photons in the control region below 1.8 GeV, with all se- 
lection criteria except the vetoes applied. Sidebands of 
the 77 mass windows are used to estimate the numbers 



of low-energy background photons that result in masses 
inside the windows. It is found that there are more such 
low-energy photons in the data than in the simulation (as 
much as 8% more at the lowest energies, below 80 MeV, 
decreasing monotonically with photon energy to approx- 
imately 2.5% above 250 MeV). 

Monte Carlo studies are used to correct for the effects 
of these differences on event-selection efficiency when 
the vetoes are imposed: —0.4% for signal events, and 
—0.9% for non-signal generic BB events. Uncertainties 
are taken to be half of the corrections. Differences be- 
tween 7r° and rj line shapes in data and simulated events 
could also potentially affect the BB efficiencies, but such 
differences proved to be negligible. 



3. Systematics of the Lepton Tag Efficiency 

There are two contributions to signal efficiency sys- 
tematics from the lepton tagging. The first is the un- 
certainty in the semileptonic branching fraction (for the 
non-signal B in the event), averaged over the lepton ac- 
ceptance for the current analysis. This is addressed in 
Sec. IVIGI and resuhs in a systematic correction and un- 
certainty of 1.047 ± 0.013. 

The second contribution arises from possible differ- 
ences between data and MC samples in the lepton iden- 
tification efficiencies. These identification efficiencies in 
the simulation are calibrated as a function of lepton mo- 
mentum to those in data using control samples of low- 
multiplicity (Bhabha and /Lt/X7) events. To measure the 
additional effect of the high-multiplicity environment in 
signal events, fitted J/ijj yields are compared in data and 
MC samples of reconstructed B J/ip {J/ip 
events, both with and without particle identification re- 
quirements applied to the leptons. This is done sepa- 
rately for e+e~ and /i^/i^ decays. The resulting system- 
atic correction factor for a single lepton, averaged over 
the mix of electron and muon tags in the current analy- 
sis, is 0.989 ± 0.004. This result is also included for the 
BB background systematics in Sec. I VI HI 

4- Other Uncertainties in Event-selection Efficiency 

A systematic uncertainty is assigned to the MC com- 
putation of the efficiency of the neural-network selection 
criteria. The control samples used to compare data and 
MC efficiencies are inclusive 7r° samples, created by ap- 
plying the standard event selection criteria to data and 
to BB background events, but with the 7r° veto inverted, 
i.e., an event is accepted if it has a 77 mass combination 
inside the veto window. The 77 mass spectra confirm 
that most of these events are due to actual 7r° produc- 
tion. Off-resonance-subtracted data are compared to the 
simulated BB sample. The efficiencies of the neural- 
network criteria for signal MC and BB background MC 
events show similar increases with E* To validate use 
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TABLE VIII: Systematic correction factors and uncertainties 
on the signal efficiency in _B — >■ Xs'y branching fraction mea- 
surements. Corrections are relative to the signal Monte Carlo 
simulation. "HE7" stands for the high-energy photon. 



Effect 




Value 


IIE7 detection efficiency 


0.9885 


± 


0.0065 


IIE7 energy scale 


1.0 


± 


0.0025 


HE7 resolution 


1.0 


± 


0.001 


HE7 lateral moment requirement 


1.0 


± 


0.003 


HE7 isolation requirement 


1.0 


± 


0.020 


Tv'^ and rj vetoes 


0.996 


± 


0.002 


Lepton PID 


0.989 


± 


0.004 


B semileptonic BF 


1.047 


± 


0.013 


Neural network 


1.0 


± 


0.012 


Hadronization model 


1.0 


± 


0.011 


Combined 


1.019 


± 


0.030 



of the 7r° control sample, neural-network output distri- 
butions for signal and control samples were compared in 
a narrow range of 1.8 < -E* < 2.0 GeV and found to be 
quite similar. Data-MC efficiency comparisons for the 
control samples are made separately for the electron and 
muon neural networks, and differences are weighted by 
the fractions of electron and muon tags in the standard 
event selection. This average difference of 1.2% is taken 
as a systematic uncertainty. 

Lastly, the signal efficiency has some small variation 
with the specific final hadronic state. The overall ef- 
ficiency is thus sensitive to whether the JETSET model 
implemented in the simulation properly describes the 
hadronization process. Measured data-MC differences 
from the BaBAR sum-of-exclusives B — ^ Xg^ analysis [4lj 
are used to reweight the hadronic multiplicity distribu- 
tion of the simulated X^ final state, and, separately, the 
fraction of final states which contain at least one 7r°. Each 
efficiency change is taken as a systematic uncertainty. 
Combining the two effects in quadrature, the total sys- 
tematic uncertainty due to modeling of the hadronization 
process is 1.1%. 



5. Overall Efficiency Systematics 

Table IVIIII summarizes the efficiency corrections and 
their estimated uncertainties. Nearly all of these effects 
are independent of photon energy E*, so the tabulated 
values apply both to wide bins and 100-MeV bins. The 
only exceptions are the small energy-scale and resolution 
uncertainties, which are folded into the yield spectrum 
(Fig. [TT|l : the Table presents the values for a bin from 
1.8 to 2.8 GeV. The correction factors are included in all 
values of efficiency quoted subsequently in this paper. 



TABLE IX: Summary of relative systematic uncertainties on 
the signal branching fraction. In addition to the contributions 
from the three factors in Eq. ([9]) (the systematic uncertainty 
on signal yield is due to that on BB background), there is 
a cross-term arising from correlations between background- 
yield and signal-efficiency uncertainties. 



E; Range (GeV) 


1.8 to 2.8 


1.9 to 2.8 


2.0 to 2.8 


Signal efficiency 


0.031 


0.031 


0.031 


BB background 


0.078 


0.051 


0.032 


Cross-terms 


0.029 


0.024 


0.019 


Count of T(45) events 


0.011 


0.011 


0.011 


Total (quadrature sum) 


0.090 


0.065 


0.050 



6. Combining Yield and Efficiency Uncertainties 

Table IIXI summarizes all systematic uncertainties for 
the branching fraction measurement. 

The fractional branching fraction uncertainty due to 
BB background is energy- dependent primarily because 
the ratio B/S of background yield to signal yield de- 
creases sharply with increasing E*. 

Similar contributions to efficiency affect the MC com- 
putations of both the BB background yield B and the 
signal efficiency Cgig , so some systematic uncertainties are 
common to both and hence are treated as correlated in 
evaluating Eq. ^ . Because of the direct calibration of tt" 
and 77 contributions to the BB background yield against 
data, some correlated effects are reduced to an insignif- 
icant level. We consider these remaining correlated ef- 
fects: 

• the systematic uncertainty due to high-energy pho- 
ton efficiency, which enters the BB yield predomi- 
nantly via the low-energy photon efficiency correc- 
tion to the tt'^ and ij components; 

• the TT^{ri) veto efficiency, which affects all back- 
ground components; 

• the semileptonic branching fraction, which for BB 
backgrounds affects lepton tags for non-Tr^/yy com- 
ponents, and also those events in which an electron 
fakes the high-energy photon signature. 

These correlated effects result in a cros s-ter m between the 
uncertainties in egig and B of 0.0178-\/B/S for the energy 
ranges considered in Table IIXI Like the BB yield con- 
tribution itself, this decreases with increasing E*. For 
the 100-MeV bins used in the spectrum measurements 
(Sec. IXII) . an additional energy dependence is allowed 
for in the semileptonic branching fraction cross-term. It 
arises because the variation of the uncertainty with lep- 
ton energy given in Fig. [TU] directly applies to the elec- 
tron backgrounds. For each of the three cross-term con- 
tributions, the uncertainty is treated as fully correlated 
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between energy bins, and the three corresponding error 
matrices are then summed. 

The total systematic uncertainty on the branching 
fraction, also given in Table ITXl is the sum in quadrature 
of the contributions from yield, efficiency, cross-terms 
and Nn-n. 



B. Model-dependence Uncertainties of the Signal 
Efficiency 

The signal efficiency egig is estimated with MC simu- 
lated spectra. The central value depends on the B — 
Xs^ model chosen, and thus has an associated model- 
dependent uncertainty. HFAG [1^ has provided world 
average values of the HQET parameters mi, and /i^ (and 
others) in the kinetic scheme, obtained from a combined 
fits to measurements of i? Xcii^ moments and previ- 
ous measurements of i? — >■ Xsj moments. (The small 
samples of earlier BaBAR B — )■ data used by HFAG 
do not lead to significant correlations between the fit re- 
sults and the data presented here.) The central values 
of the efficiency in each of three energy ranges for the 
current analysis are determined by computing efficien- 
cies for several kinetic-scheme models with to;, and /i^ 
close to the values found in the global HFAG fit, and 
interpolating to the HFAG values. For an energy range 
1.8 < E* < 2.8 GeV, the corresponding signal efficiency 
is 0.02573. 

Three considerations enter the estimate of model de- 
pendence. First, the error ellipse associated with the 
HFAG fit is used to estimate an efficiency uncertainty. 
Second, the central values for mf, and /ij from the HFAG 
fit to i? — 7> Xcii' moments only, and from a similar fit ^] 
using B — > Xs"f moments only but constraining other 
HQET parameters based on the combined fit, are con- 
sidered. The largest efficiency deviation is that from the 
B — A's7-only fit, so that is assigned as the kinetic- 
scheme uncertainty. Third, a procedure which translates 
HQET parameters from the kinetic scheme to the shape 
function scheme is applied to the combined-fit results to 
provided mb and /x^ values [25|, H^l for an efficiency es- 
timate in the shape function scheme. The difference be- 
tween that estimate and the central value in the kinetic 
scheme is added to the kinetic-scheme uncertainty (lin- 
early because both effects are systematic shifts rather 
than random variations), and taken as a symmetric un- 
certainty. Lastly, this is combined in quadrature with an 
uncertainty due to the choice of scale factor in the scheme 
translation. For the range 1.8 to 2.8 GeV, the three effects 
together yield Ae^ig = (0.00025 -|- 0.00019) © 0.00024 = 
0.00051. 

The signal efficiency, and associated model errors, for 
three photon energy ranges is given in Table [Xj 



TABLE X: 


Signal efficiency 


central values and model- 


dependence errors, for various 


ranges of measured E*. 


e; 


Range (GeV) 


^sig(%) 




1.8 to 2.8 


2.573 ±0.051 




1.9 to 2.8 


2.603 ± 0.038 




2.0 to 2.8 


2.641 ±0.029 



C. Branching Fraction Results 

Table lXll shows the branching fractions B{B — > Xs-^-dJ) 
for three ranges of measured E*, from applying Eq. © 
with the efficiencies obtained in Sec. IX Bl 

In order to compare directly to theoretical predictions, 
the measurement for each energy range in the CM frame 
is converted to a branching fraction in the corresponding 
range of true energy in the B frame. The factor a needed 
to accomplish this is determined from MC simulation us- 
ing the same methods for choosing a central value (based 
on the HFAG world-average HQET parameters in the 
kinetic scheme) and for estimating model dependence as 
are used for esig fSec. lXB]) . Values of a and the resulting 
values of B{B Xg+dl) are also presented in Table IXll 
The model-dependence uncertainties on a and 1 /csig are 
positively correlated; models with a larger fraction of the 
spectrum at low energy have larger average l/csig and 
usually larger a. Hence the fractional model-dependence 
errors on ggig and a are linearly added. Should the HFAG 
values for the kinetic-scheme parameters change in the fu- 
ture, the Appendix provides a prescription for adjusting 
both esig and a, and hence the central branching fraction 
values, for such a change. 

Finally, the contribution of B{B Xd'y) is ac- 
counted for by multiplying B{B Xs+dl) by 1/(1 + 
{\Vtd\/\Vts\)'^) = 0.958 ±0.003. This leads to the results, 
also presented in Table IXll for B{B — >■ Xgj) in true-i?-^ 
ranges, with the small additional uncertainty from this 
factor included in the systematic error. 

Because the events in the three energy ranges are 
mostly in common, even the statistical uncertainties on 
the three branching fractions are highly correlated. The 
overall correlation matrix for all statistical and system- 
atic effects (including yield-efficiency cross-terms, but ex- 
cluding model dependence) are given in Table IXlIl 

XI. UNFOLDED SPECTRUM 

The theoretical predictions of the photon energy spec- 
trum are made in the B-mcson rest frame in terms of the 
true photon energy E!^ . However the measured spec- 
trum in Fig. [11] is measured in the T(4S') frame in terms 
of the reconstructed E* after the event selection require- 
ments. To convert the measured spectrum to one that 
can be directly compared to predictions requires correct- 
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TABLE XI: Branching fractions in several photon energy ranges of both measured E* (CM frame) and true E-y {B rest frame), 
along with the adjustment factor a between them. Uncertainties on branching fractions are statistical, systematic and model- 
dependence, respectively. The error on the adjustment factor a is a model-dependence uncertainty, treated as fully correlated 
with that on the initial B. 



Energy Range 


BiB ^ X,+d7) (10-*) 
in Measured E* Range 


Factor a 
to True E^ 


B{B ^ X,+d7) (10-*) 
in True E^ Range 


BiB X,j) (10-*) 
in True E-, Range 


1.8 to 2.8 GeV 

1.9 to 2.8 GeV 
2.0 to 2.8 GeV 


3.271 ± 0.154 ± 0.294 ± 0.065 
3.019 ± 0.136 ± 0.196 ± 0.044 
2.745 ± 0.120 ± 0.137 ± 0.030 


1.0233 ± 0.0042 
1.0356 ± 0.0045 
1.0657 ± 0.0045 


3.347 ± 0.158 ± 0.301 ± 0.080 
3.126 ± 0.141 ± 0.203 ± 0.059 
2.925 ± 0.128 ± 0.146 ± 0.045 


3.207 ± 0.151 ± 0.288 ± 0.077 
2.995 ± 0.135 ± 0.194 ± 0.057 
2.802 ± 0.122 ± 0.140 ± 0.043 



TABLE XIL The correlation matrix for the measured branching fractions in three energy ranges, including all statistical and 
systematic but not model-dependence uncertainties. 



E-f Range 


1.8 to 2.8 GeV 


1.9 to 2.8 GeV 


2.0 to 2.8 GeV 


1.8 to 2.8 GeV 


1.00 


0.94 


0.84 


1.9 to 2.8 GeV 




1.00 


0.92 


2.0 to 2.8 GeV 






1.00 



ing for selection efficiency and detector acceptance, and 
unfolding two resolution effects. These are detector res- 
olution and Doppler smearing of the photon energy. The 
transformation of the measured E* spectrum to an 
spectrum thus requires four steps: 

1. Correcting for the event selection efficiency. 

2. Unfolding the effects of detector resolution. 

3. Correcting for the detector acceptance. 

4. Unfolding the Doppler smearing. 

Each of these steps requires the use of the MC sim- 
ulation to either estimate the efficiency and acceptance 
or model the resolution and smearing. The effects of 
calorimeter resolution and Doppler smearing on the pho- 
ton spectrum are unfolded using a simplified version of an 
iterative method [61J. This simplified method has been 
used previously by the BaBAR collaboration in a mea- 
surement of the e+e" — )■ 7r+7r-(7) cross section [g^I- An 
introduction to this method is followed by a description 
of the implementation used here and then the results and 
systematic uncertainties. The notation used for the pho- 
ton energy is: 

• E* is the energy measured in the T(4S') rest frame 
after event selection. 

• E* is the true photon energy in the T(45') rest 
frame. Its spectrum is obtained after steps 1-3 
above. 

• E^ is the true photon energy in the _B-meson rest 
frame. Its spectrum is obtained after steps 1-4 
above. 



A. Overview of the Unfolding Technique 

The effects of detector resolution and Doppler smear- 
ing each require a separate unfolding but the procedure 
for each is identical. In this overview the unfolding of 
the detector resolution is described. The unfolding of the 
Doppler smearing uses the same procedure. First some 
general considerations for the unfolding are given before 
describing the features of implementation used. 

The spectrum is measured in twelve 100-MeV bins be- 
tween 1.6 and 2.8 GeV and one 70-MeV bin between 1.53 
and 1.6 GeV. The detector resolution can cause a migra- 
tion between bins which is described by a transfer matrix 
A, whose elements Aij are the number of events gener- 
ated in bin j that are reconstructed in bin i. Identical 
binning is used for the generated and smeared spectra 
so that Aij is a square matrix. The transfer matrix is 
derived from MC simulation using an assumed model for 
the spectrum. It is then used to construct a folding ma- 
trix Pij and an unfolding matrix 



' " A■^ ' 

where Pij is the probability for an event generated in 
the bin j to be reconstructed in the bin i and Pij is the 
probability of the reconstructed event in the bin i coming 
from the generated bin j. N is the number of bins. In 
principle the unfolding matrix can now be directly ap- 
plied to the reconstructed spectrum to unfold the resolu- 
tion effects. There are, however, two significant problems 
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with this approach. The first is that it assumes the sim- 
ulated model perfectly describes the data. The second 
is that any significant statistical fluctuations in the re- 
constructed spectrum can be unfolded into several bins, 
causing unstable and unreliable results. 

The technique adopted mitigates these problems. It 
begins by simulating an approximate model of the spec- 
trum that is normalized to data in the range 1.8 < E* < 
2.8 GeV. This model is referred to as the initial model. 
The difference between this model and the data is then di- 
vided into two parts. The first part is attributed to a gen- 
uine difference between the model and the true data spec- 
trum and is used to modify the transfer matrix, equiv- 
alent to changing the initial model. The second is at- 
tributed to statistical and systematic fluctuations and is 
not unfolded using the unfolding matrix, but rather used 
to correct the model spectrum so that significant fiuc- 
tuations in the reconstructed spectrum are propagated 
to the unfolded true data spectrum. The division of the 
difference into these two parts is accomplished using a 
bin-dependent regularization function / with a tunable 
parameter A. The value of / varies from to 1 according 
to the value of A so that a fraction / comprises the true 
model-to-data difference and a fraction 1 — / the statisti- 
cal and systematic fluctuation. A priori the value of A is 
unknown, but can be estimated using an MC technique 
described in Sec. IXIBI 

The technique has been tested extensively in simulated 
data and found to give reliable and stable results. 

B. Implementation of the Unfolding 

This nominal initial model is found by comparing the 
data to a set of models in the kinetic, shape function, 
and KN schemes using different values of HQET param- 
eters. Each model is passed through the full simulation 
and event selection. Figure [T^] shows the comparison of 
the data to a range of models that describe the data at 
the one-cr level. The closest match is chosen by construct- 
ing a function formed from the bin-by-bin differences 
of the data and the generated spectrum using 100-Mev 
bins in the signal range 1.8 < -E* < 2.8 GeV and the full 
covariance matrix. It is found that a model in the shape 
function scheme with {nib — 4.51 GeV, /i^ = 0.46 GeV^) 
best describes the data. The other models shown in 
Fig. [12] are used to optimize the A parameters for the 
two unfolding steps and to estimate mo del- dependence 
systematic uncertainties. 

The unfolding method begins by correcting the mea- 
sured data spectrum for selection efficiency in each bin. It 
is then compared with the reconstructed simulated spec- 
trum of the initial model by computing the difference 

Adi = di — Cvi . 

Here di is the number of efficiency-corrected recon- 
structed data events in the z*'' bin, r.^ is the number of 
efficiency-corrected reconstructed simulated events and 



C normalizes the initial model spectrum to the data in 
the signal range 1.8 < i?* < 2.8 GeV. A fraction / of Ad, 
comes from a true difference between the model and the 
data spectrum, while the remaining fraction 1 — / is due 
to a fluctuation in cither the signal or in the background 
subtraction. The function / is a regularization function 
with a tunable parameter A: 

/(Ad„a„A) = l-e"^^)' , 

where ai is the error in di. There are several choices of 
regularization functions suggested in reference [6l| . Each 
function has the property that it varies monotonically 
between and 1 as the combination Adi/{Xai) changes 
from to 00. The procedure is found to be insensitive 
to the particular choice, so the simplest is chosen. The 
value of the regularization parameter A thus determines 
the fraction / of the difference that is unfolded in each 
bin. 

An ensemble of 40,000 simulated model spectra is used 
both to optimize A and to derive the error matrix for the 
unfolded spectrum. These spectra have been generated 
using an error matrix that is constructed from the errors 
in Table llVi the bin-to-bin correlations in Table |V] and 
the correlations between the background subtraction and 
the efficiency systematics described in Sec. IX A 61 Each 
spectrum is unfolded with a range of values of A and 
then corrected for acceptance. The error matrix of the 
unfolded spectra, O, is computed from the ensemble us- 
ing the output distributions of energy in each bin and the 
correlations between these distributions. This is then re- 
peated using different models to construct the unfolding 
matrix. A function is formed using a vector of the 
unfolded yields, the inverse of the error matrix and 
a vector of the true value of the original generated MC 
spectra t: 

Only bins in the signal region 1.8 < E* < 2.8 GeV are 
used for the optimization. The function is then used 
to find the value of A that most closely reproduces t for 
all models. The optimal value of A is then used to unfold 
the data. The detector resolution unfolding is performed 
using A = 0.5. 

The unfolding matrix used to unfold the data P/^- is 
constructed from a modified transfer matrix A',-. It is 
modified by adding the folded difference between the ini- 
tial model and the data: 

= CA,, +AdlP,, 
p' = ^''^ 

The unfolded data spectrum Uj is then obtained from 

N 

^3 = Ci, -K^{/(Ad,;,(Td,,A)-Ad, 

1=1 

+ [l-/(Ad„ad,,A)] • Ad,-5„} . 
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It is expressed in terms of a correction to the true value 
of the initial model {Ctj). The second term is that part 
of the difference between the initial model and the data 
that is to be unfolded using the unfolding matrix P'ij. 
The third term is the part of the difference attributed to 
statistical or systematic fluctuation that is not unfolded. 
This procedure was iterated, but its output was found 
to have converged after just one application, so that first 
iteration provides the results presented below. The pro- 
cedure for unfolding the Doppler smearing is identical 
except that the optimal value of A = 1.0. 

C. Results of the Unfolding 

The measured spectrum shown Fig. [TT] and the cor- 
responding yields and uncertainties in Table IIVI are the 
starting point for the unfolding. First the spectrum is 
corrected for the selection efficiency, taking into account 
the additional correlated errors between the efficiency 
and the background estimation described in Sec. IX A 61 
Then the resolution smearing is unfolded and the resul- 
tant spectrum corrected for detector acceptance to give 
a spectrum in bins of presented in Table IXlIIl 

The estimation of the statistical, systematic and model- 
dependence uncertainties is described in Sec. IXI Dl To 
provide complete yield uncertainties, the 3.1% energy- 
independent uncertainty on efficiency (see Table lIXp is 
included in the systematic uncertainty. The spectrum is 
shown in Fig. [131 

The Doppler smearing is then unfolded starting from 
Fig. [13] and Table IXlIIl The resulting yields in bins of 
are converted to partial branching fractions by di- 
viding by the number of B mesons in the on-resonance 



TABLE XIII: The unfolded E' spectrum and its uncer- 
tainties in numbers of produced events. The total error is the 
sum in quadrature of the statistical, systematic and model- 
dependence errors 



i;;"'"" (GeV) 


Yield 




Error (events) 




Stat. 


Syst. 


Model 


Total 


1.53 


to 


1 


60 


24620 


15193 


24749 


657 


29115 


1.60 


to 


1 


70 


51556 


12190 


21593 


365 


25140 


1.70 


to 


1 


80 


4244 


10631 


15598 


42 


18427 


1.80 


to 


1 


90 


22346 


8999 


11278 


208 


14217 


1.90 


to 


2 


00 


22506 


7252 


7565 


94 


10626 


2.00 


to 


2 


10 


22177 


5705 


4708 


1512 


7461 


2.10 


to 


2 


20 


27518 


4773 


2865 


939 


5406 


2.20 


to 


2 


30 


42298 


4140 


2037 


1073 


4673 


2.30 


to 


2 


40 


39193 


4010 


1542 


1256 


4384 


2.40 


to 


2 


50 


43214 


3755 


1671 


1140 


4164 


2.50 


to 


2 


60 


29488 


3560 


1065 


611 


3789 


2.60 


to 


2 


70 


20025 


2784 


723 


75 


2857 


2.70 


to 


2 


80 


610 


2446 


179 


141 


2467 
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FIG. 13: The photon spectrum after unfolding the effects 
of calorimeter resolution and correcting for the selection effi- 
ciency and detector acceptance. The inner error is statistical 
only, the outer includes statistical, systematic and model- 
dependence errors added in quadrature. The vertical line 
shows the boundary between the lower control region and the 
signal region. 



TABLE XIV: Partial branching fraction in bins of ob- 
tained from the unfolded spectrum. These values describe 
the shape of the spectrum in the B rest frame and provide a 
crosscheck (see Sec. XI. E) of the integrated branching frac- 
tions, but are not intended as primary branching fraction re- 
sults. (The integrated branching fractions reported in Ta- 
ble IXII are more precise and less susceptible to bieis, as ex- 
plained in Sec lVIlH ') The total error is the sum in quadrature 
of the statistical, systematic and model-dependence errors. 
The model error is relatively large in the bins above 2.4 GeV, 
but anticorrelated between neighboring bins, as discussed in 
Sec. IXIDI Hence combined 200-MeV bins for this region are 
shown at the bottom of this table and in Fig. 1141 



(GeV) 




d7) 


Error 




(10-') 


Stat 


Syst 


Model 


Total 


1.50 to 1.60 


2.53 


1.59 


2.52 


0.33 


2.97 


1.60 to 1.70 


7.76 


1.95 


3.90 


0.31 


4.44 


1.70 to 1.80 


0.25 


1.53 


2.07 


0.06 


2.48 


1.80 to 1.90 


2.81 


1.30 


1.45 


0.03 


1.87 


1.90 to 2.00 


3.16 


1.05 


1.03 


0.10 


1.45 


2.00 to 2.10 


2.67 


0.83 


0.65 


0.28 


1.06 


2.10 to 2.20 


3.56 


0.70 


0.38 


0.16 


0.76 


2.20 to 2.30 


5.44 


0.60 


0.28 


0.26 


0.69 


2.30 to 2.40 


5.37 


0.58 


0.23 


0.16 


0.62 


2.40 to 2.50 


5.80 


0.53 


0.24 


0.99 


1.13 


2.50 to 2.60 


6.46 


0.59 


0.26 


0.80 


1.02 


2.60 to 2.70 


0.00 


0.11 


0.01 


0.12 


0.16 


2.70 to 2.80 


-0.12 


0.21 


0.01 


0.10 


0.23 


2.40 to 2.60 


12.25 


0.79 


0.47 


0.19 


0.92 


2.60 to 2.80 


-0.12 


0.24 


0.01 


0.22 


0.32 
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TABLE XV: The correlation matrix for the errors on signal yields shown in Table IXlIIl in bins of E* Statistical, systematic, 
and model errors are included. Columns are labeled by the value of E* at the lower edge of the bin. 



e; (Gcv) 


1.53 


1.6 


1.7 


1.8 


1.9 


2.0 


2.1 


2.2 


2.3 


2.4 


2.5 


2.6 


2.7 


1.53 to 1.6 


1.00 


0.65 


0.63 


0.60 


0.55 


0.44 


0.29 


0.16 


0.11 


0.13 


0.07 


0.03 


0.04 


1.6 to 1.7 




1.00 


0.63 


0.59 


0.55 


0.46 


0.31 


0.20 


0.09 


0.15 


0.06 


0.09 


0.05 


1.7 to 1.8 






1.00 


0.57 


0.52 


0.41 


0.29 


0.16 


0.10 


0.12 


0.04 


0.04 


0.03 


1.8 to 1.9 








1.00 


0.50 


0.39 


0.26 


0.15 


0.09 


0.17 


0.05 


0.09 


0.01 


1.9 to 2.0 










1.00 


0.38 


0.25 


0.14 


0.09 


0.12 


0.05 


0.01 


0.05 


2.0 to 2.1 












1.00 


0.25 


0.18 


0.02 


0.04 


0.03 


0.03 


-0.02 


2.1 to 2.2 














1.00 


0.16 


0.00 


0.07 


0.06 


0.03 


0.01 


2.2 to 2.3 
















1.00 


0.00 


0.03 


0.06 


0.08 


-0.02 


2.3 to 2.4 


















1.00 


0.17 


0.07 


0.04 


0.05 


2.4 to 2.5 




















1.00 


0.11 


0.05 


0.05 


2.5 to 2.6 






















1.00 


0.02 


-0.06 


2.6 to 2.7 
























1.00 


0.02 


2.7 to 2.8 


























1.00 



TABLE XVL The correlation matrix for the errors on partial branching fractions shown in Table IXlVl in bins of E^^ . Statistical, 
systematic, and model errors are included. Columns are labeled by the value of E!^ at the lower edge of the bin. 



E^ (GeV) 


1.53 


1.6 


1.7 


1.8 


1.9 


2.0 


2.1 


2.2 


2.3 


2.4 


2.6 


1.53 to 1.6 


1.00 


0.51 


0.53 


0.53 


0.50 


0.36 


0.20 


0.05 


0.07 


0.13 


0.18 


1.6 to 1.7 




1.00 


0.61 


0.57 


0.55 


0.48 


0.28 


0.15 


0.01 


0.12 


0.08 


1.7 to 1.8 






1.00 


0.54 


0.48 


0.41 


0.24 


0.11 


0.03 


0.07 


0.10 


1.8 to 1.9 








1.00 


0.48 


0.35 


0.21 


0.08 


0.01 


0.15 


0.11 


1.9 to 2.0 










1.00 


0.37 


0.15 


0.05 


0.05 


0.11 


0.16 


2.0 to 2.1 












1.00 


0.27 


0.12 


-0.07 


0.04 


-0.12 


2.1 to 2.2 














1.00 


0.19 


-0.14 


0.07 


-0.09 


2.2 to 2.3 
















1.00 


0.00 


-0.03 


-0.23 


2.3 to 2.4 


















1.00 


0.11 


0.20 


2.4 to 2.6 




















1.00 


0.16 


2.6 to 2.8 






















1.00 



data sample, 2Ngj^. These branching fractions are pre- 
sented in Table IXTVl An additional 1.1% has been in- 
cluded in the systematic error to account for the uncer- 
tainty in Ngg. Figure [T4l shows this photon spectrum in 
the B rest frame. The spectrum is compared to that for 
a kinetic-scheme model with parameters = 4.60 GeV 
and /ij — 0.45 GeV^, close to HFAG world averages. 
(The if* (892) has not been substituted for the highest- 
energy part of the spectrum, because the unfolded data 
cannot resolve such a peak.) 

The correlation matrices corresponding to Table IXIIII 
and IXIVI are given in Tables |XV] and IXVll respectively. 
These matrices have a complex structure because many 
effects contribute. At low energies (the upper left quad- 
rant) they are dominated by the highly-correlated uncer- 
tainties in the BB backgrounds. At higher energies, the 
uncorrelated statistical uncertainty is relatively more im- 
portant, along with the smaller but fully-correlated sys- 



tematic uncertainty on efficiency. The contributions from 
the unfolding itself and from model dependence can be 
negative. Hence in the lower-right quadrant, where other 
correlations are weak, the net result can be close to zero 
or negative. 

The numbers in Tables IXI VI and IX VII can be used to fit 
the measured spectral shape to any theoretical prediction 
in the _B-meson rest frame. 



D. Statistical, Systematic and Model-Dependence 
Uncertainties in the Unfolding 

The dominant uncertainty in the bins of the unfolded 
spectrum is due to the BB subtraction described in 
Sec. IVIII The statistical and systematic errors on the 
efficiency-corrected yields are propagated using the en- 
semble MC technique described previously. A num- 
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FIG. 14; The photon spectrum after unfolding resolution and 
Doppler smearing, shown as a partial branching fraction (AB, 
see Table lXTVl caption'). The inner error is statistical only, the 
outer includes statistical, systematic and model-dependence 
errors added in quadrature. Section fXI DI explains why results 
above 2.4 GeV are shown in wider bins. The vertical line 
shows the boundary between the lower control region and the 
signal region. The curve shows the spectrum in a kinetic- 
scheme model (see text), normalized to the data from 1.8 to 
2.8 GeV/c. 



TABLE XVII: The change in the number of events in each 
bin of the unfolded photon spectrum after shifting the photon 
energy scale by ±0.3%. The absolute value of the largest 
difference (4- or — ) is shown after resolution unfolding {E^ '"^"^ 
bins) and both resolution and Doppler smearing unfolding 
{E^ bins). In both cases efficiency and acceptance corrections 
have been applied. These changes are included in the final 
systematic errors in Tables IXIIII and IXIVI assuming 100% 
correlation between the bins. 



smearing the calorimeter resolution in the MC simula- 
tion by an additional 1%, as determined by data com- 
parisons in Sec. IX A II The only significant effects are 
found to be in the photon energy scale shift. Table IXVlII 
shows the bin-by-bin change in the event yields due to the 
photon energy shift. For each bin, the absolute value of 
the largest difference (-1- or — ) is taken as the systematic 
uncertainty, and 100% bin-to-bin correlation is assumed. 
This error is combined in quadrature with the systematic 
error propagated from the measured E* spectrum and is 
included in Tables Km] and |XlVl 

To assess the model dependence, the unfolding is per- 
formed with a range of models. In each case the same 
model is used for the entire procedure including effi- 
ciency and acceptance corrections, and unfolding the de- 
tector resolution and Doppler smearing. Figure [T^ shows 
two models that could plausibly describe the data at 
the one-sigma level. These are a shape function model 
with (mb = 4.40 GeV, /i^ = 0.52 GeV^) and a kinetic- 
scheme model {nib = 4.60 GeV, = 0.45 GeV^,^^ = 
0.27 GeV^). To set the model-dependence error we un- 
fold the nominal simulated model (shape function: nib — 
4.51 GeV, fil = 0.46 GeV^) with one of these two models. 
The largest bin-by-bin difference is taken as the model 
error in the unfolded spectrum with 100% correlation be- 
tween each bin. The model-dependence error is generally 
much smaller than the systematic error except for the 
unfolding of the Doppler smearing close to the kinematic 
limit {E^ w mB/2). The steeply falling spectrum at this 
limit leads to a much greater sensitivity to the model, 
which results in a large error that is anti-correlated be- 
tween the 2.4-2.5 GeV and 2.5-2.6 GeV bins. To avoid 
this edge effect the two bins are summed. This is also 
done for the 2.6-2.8 GeV range. 



Energy Range (GeV) 


Change (events) 


E; "^"^ Bins 


e!^ Bins 


1.53 to 1.60 


222.1 


220.2 


1.60 to 1.70 


190.6 


191.0 


1.70 to 1.80 


261.1 


261.6 


1.80 to 1.90 


354.4 


354.8 


1.90 to 2.00 


493.2 


492.0 


2.00 to 2.10 


622.9 


622.2 


2.10 to 2.20 


640.3 


658.5 


2.20 to 2.30 


428.4 


461.1 


2.30 to 2.40 


528.7 


598.9 


2.40 to 2.50 


1184.2 


1292.5 


2.50 to 2.60 


1080.6 


967.6 


2.60 to 2.70 


490.8 


475.7 



E. Crosscheck of Branching Fraction 



ber of possible uncertainties in the unfolding procedure 
were considered. These included changing the regular- 
ization parameter A to zero which changes / to 1.0 in 
all bins, changing the normalization factor C accord- 
ing to the 10% uncertainty in the measured value of 
B{B — > Xsj), varying the energy scale by ±0.3%, and 



The numbers in Table IXIVI are used to obtain inte- 
grated branching fractions B{B Xsj) for purposes 
of comparison with the reported results from Sec. IX Gl 
The AB{B — > Xg+dj) values are summed over E^^ inter- 
vals, with the errors combined including correlations (Ta- 
ble |XVT|. Lastly, a factor of 0.958 is applied to account 
for the B — > X^'j contribution. As explained in Sec lVIIll 
the unfolded yields are based on a different choice of 
model than that used to extract the B{B — )■ Xgj) re- 
sults for this analysis, and hence are not intended to be 
used for such results. This procedure has been carried 
out for one energy range, 1.8 < E^^ < 2.8 GeV. 

There are two contributions to the uncertainty beyond 
those implied by Tables IXIVI and IXVII First, there is 
the small (1.1%) uncertainty on Ng^. Second, because 
the range of models used to estimate mo del- dependence 
uncertainty is data-driven, that uncertainty is positively 
correlated with the systematic uncertainty on the signal 
yield. This gives a B{B Xsj) for 1.8 < E^ < 2.8 GeV 
of (3.36±0.19±0.34±0.08) x lO"" = (3.36±0.43) x 10-^ 
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where the first set of errors are statistical, systematic and 
model, and their combination in the second form takes 
the model-systematic correlation into account. 

This value may be compared to the reported branch- 
ing fraction of (3.20 ± 0.15 ± 0.29 ± 0.08) x 10""* = 
(3.20 ± 0.33) X 10"-* from TableKU the three uncertain- 
ties (independent in that case) are added in quadrature. 
The difference in the central values is due to the differ- 
ent choice of the central model; if a data-like model had 
been used in Sec. [XJ the extracted branching fraction 
would have been 3.36 x 10~^, the same value obtained 
with unfolding. 

The smaller statistical and systematic uncertainties on 
the branching fraction from Table IXII are in large part a 
consequence of applying the efficiency correction to a sin- 
gle wide bin of photon energy. As discussed in Sec. IVIIIl 
this de-emphasizes the importance of the uncertainties in 
the lowest-energy region, where signal efficiency is lowest 
and background uncertainties are largest. The branching 
fraction as derived from the unfolded spectrum of neces- 
sity relies upon efficiency corrections in 100-MeV bins. In 
addition, the combined uncertainty on the latter result is 
increased by the model-background correlation discussed 
above, an effect which does not occur when the model 
range is chosen as described in Sec. IX Bl 

F. Moments of the Spectrum 

The moments of the spectrum provide information to 
measure the HQET parameters mt and /xj in the kinetic 
scheme [l^. The first, second and third spectral mo- 
ments, Ei,E2,E3 are defined in Sec. HI Eq. ([T]). They are 
measured for three photon energy ranges: 1.8 to 2.8 GeV, 
1.9 to 2.8 GeV and 2.0 to 2.8 GeV. The moments are com- 
puted directly from the unfolded spectrum in 100-Mev 
bins given in Tables IXIIII and IXIVI using the correlation 
matrices given in Tables IXVI and IXVII 

The behavior of the moments for different photon en- 
ergy ranges has been studied theoretically in the kinetic 
scheme. The spectral moments in E* are given in 
Table IXVIIII The correlations between the moments are 
given in Table IXXI to allow fits to predictions of the mo- 
ments. The E:^ spectral moments and correlations be- 
tween the moments are given in Tables IXIXI and IXXII 

XII. CONCLUSIONS 

In summary, the B — > Xg+dJ photon energy spectrum 
in the CM frame has been measured in 347.1 fb~^ of data 
taken with the BaBAR experiment. It is used to extract 
measurements of the direct CP asymmetry for the sum of 
B — >■ Xs'y and B — > X^j, the branching fraction for B — > 
Xs7, and the spectral shape and its energy moments in 
the _B-mcson rest frame. The result for CP asymmetry 
is 

Acp = 0.057 ± 0.060(stat) ± 0.018(syst) . 



\ \ 1 r 

SM 
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Babar hadron Tag -A- 



CLEO lepton Tag ■ ^ 



This Analysis 



-0.4-0.3-0.2-0.1 0.1 0.2 0.3 0.4 

Acp(B^X,,,y) 

FIG. 15: Measurements of Acp {B — >■ Xs+d^), with statistical 
and systematic errors. The three published results, top to 
bottom, are from references [12], [13] and [s^], respectively. 



The branching fraction and moments are presented for 
three ranges of the photon energy in the _B-meson rest 
frame, 1.8, 1.9 and 2.0 to 2.8 GeV (Tables E] and IXTX]) . 
For example, in the 1.8 to 2.8 GeV range: 

B{B X,-f) = (3.21 ± 0.15 ± 0.29 ± 0.08) x 10"^ , 
El = (2.267 ± 0.019 ± 0.032 ± 0.003) GeV and 
E2 ^ (0.0484 ± 0.0053 ± 0.0077 ± 0.0005) GeV^ , 

where the errors are from ststistics, systematics and 
model dependence, respectively, and the moments are 
defined in Eq. [T] 

Figure [T51 compares the measured Acp {B — > Xg+dl) 
to previous measurements and to the SM prediction. No 
asymmetry is observed, consistent with SM expectation. 
The current measurement is the most precise to date. 

Figure compares the measured branching fraction 
to previous measurements performed for different Ej 
ranges. This measurement supersedes the previous fully 
inclusive measurement from BaBAR. It is consistent with 
previous measurements and of comparable precision to 
the recent Belle measurement [40i. In order to com- 
pare with theoretical predictions the measurement for 
E^ > 1.8 GeV can be extrapolated down to 1.6 GeV using 
a factor provided by the HFAG collaboration 25]. They 
fit results from previous measurements of -B — > ^"^7 and 
B — )■ Xciiy to predictions in the kinetic scheme to yield 
average values of mi, and /i^. These parameters are then 
used to generate a S — > model in the kinetic scheme 
which gives an extrapolation factor of 1/(0.968 ± 0.006). 
When applied to the present result this gives B{B ^■ 
Xs-/) = (3.31±0.16±0.30±0.10)xl0-4(£;^ > 1.6GeV) 
which is in excellent agreement with the SM prediction 
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FIG. 16: The measured branching fraction for this anal- 
ysis (•) compared to previous measurements for different 
Ej ranges (minimum energies given on the left axis). 
The prev ious measurements are from CLEO (*) ^33], Belle 
(A) J40|, and BABAR using the semi-inclusive technique 
(□) [4l[. Error bars show total uncertainties. 



The various theoretical uncertainties are assumed to have 
Gaussian distributions, and are combined in quadrature. 
A point is then excluded if the negative la deviation of 
the prediction hes above the 95% confidence-level upper 
limit of the measured branching fraction extrapolated to 
1.6 GeV. The region M^± < 327 GeV is excluded at the 
95% confidence level, independent of tan /3. 



2.0GeV 



1.9GeV 



I.SGeV 



1.7GeV 




2.3 2.35 

<E,^(GeV) 



2.4 



2.0GeV 



Bm ^ X,j) = 3.15 ± 0.23 X 10~'^{E^ > 1.6 GeV) 
[15| and can be used to provide stringent constraints on 
new physics. An example is shown in Figure 1171 The 
effects of a type-II two-Higgs-doublet model (THDM) 
on B{B — >• X^j) at next-to-leading order are presented 
in Refs. [l^ |6^. Software provided by the author of 
Ref. [i^ computes an excluded region, following a pro- 
cedure described in Ref. (63 |. The branching fraction, 
including both the SM and the THDM contributions, is 
calculated for each point in the Mh± vs. tan/3 plane. 
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FIG. 17: The shaded area shows the excluded region (at the 
95% confidence level) in charged Higgs mass vs. tan /? for a 
type-II two-Higgs-doublet model, using the measured value 
of BiB Xs-f) = (3.31 ± 0.16 ± 0.30 ± 0.10) x 10"'' {E^ > 
1.6 GeV) from this analysis. This plot is based on predictions 
in references [isj and 
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FIG. 18: The measured first (top) and second (bottom) mo- 
ments from this analysis (•) compared with the previous mea- 
surement for different E^ ranges (minimum energies given 
on the left axis). These previous measurements are CLEO 
(*) [il, BABAR semi-inclusive (□) [13], and Belle (A) 
Error bars show total uncertainties. 



The effects of detector resolution and Doppler smear- 
ing are unfolded to present the photon spectrum in the 
B-meson rest frame for the first time in Fig. [M] This 
spectrum may be used to extract information on HQET 
parameters in two ways. First, the full covariance matrix 
is provided to allow any theoretical model to be fit to 
the entire spectrum. Secondly the moments have been 
extracted and can be compared to predictions for differ- 
ence energy ranges. Figure [18] compares the measured 
moments to previous measurements. 
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TABLE XVIII: The E*^''"' spectral moments and errors (± statistical ± systematic ± model-dependence). Moments are 
defined by Eq. ([ij in Sec. H 



E* Range (GeV) 


El (GeV) 


E2 (GeV2) 


Es ( GeV^) 


1.8 to 2.8 


2.275 ± 0.018 ± 0.032 ± 0.003 


0.0546 ± 0.0049 ± 0.0074 ± 0.0005 


-0.0031 ± 0.0011 ± 0.0013 ± 0.0004 


1.9 to 2.8 


2.314 ± 0.013 ± 0.017 ± 0.004 


0.0417 ± 0.0032 ± 0.0028 ± 0.0003 


-0.0013 ± 0.0007 ± 0.0005 ± 0.0003 


2.0 to 2.8 


2.350 ± 0.010 ± 0.008 ± 0.005 


0.0317 ± 0.0022 ± 0.0010 ± 0.0005 


0.0001 ± 0.0005 ± 0.0002 ± 0.0002 



TABLE XIX: The E!^ spectral moments and errors (± statistical ± systematic ± model-dependence). Moments are defined 
by Eq. ^ in Sec. H] 



Range (GeV) 


El (GeV) 


E2 (GeV^) 


E3 (GeV^) 


1.8 to 2.8 


2.267 ± 0.019 ± 0.032 ± 0.003 


0.0484 ± 0.0053 ± 0.0077 ± 0.0005 


-0.0048 ± 0.0011 ± 0.0011 ± 0.0004 


1.9 to 2.8 


2.304 ± 0.014 ± 0.017 ± 0.004 


0.0362 ± 0.0033 ± 0.0033 ± 0.0005 


-0.0029 ± 0.0007 ± 0.0004 ± 0.0002 


2.0 to 2.8 


2.342 ± 0.010 ± 0.008 ± 0.005 


0.0251 ± 0.0021 ± 0.0013 ± 0.0009 


-0.0013 ± 0.0005 ± 0.0002 ± 0.0001 



TABLE XX: The correlation matrix of the E* spectral moments. Superscripts denote the lower end of the energy range 
in GeV. 





El-^ 


El-^ 


El-' 


El-' 


El-' 


El-' 


El" 


El-' 


Ef 




1.00 


-0.88 


-0.09 


0.84 


-0.68 


-0.26 


0.61 


-0.29 


-0.30 






1.00 


-0.27 


-0.58 


0.71 


0.19 


-0.25 


0.43 


0.24 


Ek' 






1.00 


-0.29 


0.19 


0.55 


-0.23 


0.26 


0.35 


El-' 








1.00 


-0.75 


-0.31 


0.75 


-0.28 


-0.34 


El-' 










1.00 


0.12 


-0.17 


0.64 


0.30 














1.00 


-0.61 


0.62 


0.85 
















1.00 


0.25 


-0.41 


El-° 
















1.00 


-0.50 


Er 


















1.00 



TABLE XXI: The correlation matrix of the E!^ spectral moments. Superscripts denote the lower end of the energy range in 
GeV. 





El-' 


El' 


El-' 


El' 


El' 


El' 


El" 


El" 


El" 




1.00 


-0.90 


0.10 


0.84 


-0.73 


-0.05 


0.53 


-0.46 


-0.09 


El-^ 




1.00 


-0.35 


-0.60 


0.73 


-0.07 


-0.21 


0.45 


0.12 


El-' 






1.00 


-0.23 


0.14 


0.48 


-0.29 


0.30 


0.26 


El' 








1.00 


-0.82 


-0.05 


0.68 


-0.50 


-0.08 


El' 










1.00 


-0.11 


-0.27 


0.59 


0.14 


El-' 












1.00 


-0.50 


0.52 


0.59 
















1.00 


-0.59 


-0.06 


El" 
















1.00 


0.20 


El" 


















1.00 
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Appendix A: Parameterization of Branching 
Fraction Factors 

The central values of the partial branching fractions re- 
ported in Section [X CI depend on the signal efficiency ggig 



and adjustment factor a computed for a kinetic-scheme 
model with parameters nib and /i^ set to HFAG world- 
average values. This Appendix provides functional forms 
for the dependence of egig and a on these HQET parame- 
ters. In the event of possible future changes in the HFAG 
values, the information presented here would allow for a 
corresponding adjustment of the branching fraction cen- 
tral values. Since each partial branching fraction in the 
B rest frame is proportional to a/egig, the adjustment 
would be made by dividing out that combination com- 
puted for the current values of HQET parameters and 
multiplying by the same combination computed for the 
new values. 

These functions have no physical significance. They 
result from fits to the egig and a values computed by MC 
simulation for a wide range of parameters. For a grid of 
models spanning 4.5 < nib < 4.7GeV/c^ and 0.3 < n"^ < 
0.7GeV/c, these fits have fractional accuracy of better 
than 0.2% for egig and 0.1% for a. The functional form 
used for both quantities is 

f{mb,fil) = fo 

+ Mnib - 4.6 GeV/c^) + ^(v^ - 0.6 GeV/c) 

+ fsimb - 4.6 GeV/c2)( - 0.6 GeV/c) 
+ /4(mfc-4.6GeV/c2)2 , (Al) 
where the coefficients /g through have appropriate 
units to make each term diniensionless. Table IXXIII gives 
the values of these coefficients for egjg and a for each of 
the three photon-energy ranges in which branching frac- 
tions are reported. 
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